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ABSTRACT 



< 

Si 

O ' We identify satellites of isolated galaxies in the Sloan Digital Sky Survey and ex- 

' amine their angular distribution. Using mock catalogues generated from cosmological 

./V-body simulations, we demonstrate that the selection criteria used to select isolated 
galaxies and their satellites in large galaxy redshift surveys must be very strict in order 
to correctly identify systems in which the primary galaxy dominates its environment. 
£SJ ■ We demonstrate that the criteria used in many previous studies instead select pre- 

' dominantly group members. We refine a set of selection criteria for which the group 

\ contamination is estimated to be less than 7% and present a catalogue of the resulting 

. sample. 

The angular distribution of satellites about their host is biased towards the major 
axes for spheroidal galaxies and probably also for red disc galaxies (the "intermediate" 
\Q • class of Bailin & Harris), but is isotropic for blue disc galaxies, i.e. it is the colour of the 

C> host that determines the distribution of its satellites rather than its morphology. The 

, similar anisotropy measured in this study as in studies that were dominated by groups 

implies that group-specific processes are not responsible for the angular distribution. 
Satellites that are most likely to have been recently accreted, late-type galaxies at 
large projected radii, show a tendancy to lie along the same axis as the surrounding 
/\ ' large scale structure. The orientations of isolated early and intermediate-type galaxies 

also align with the surrounding large scale structures. 

We discuss the origin of the anisotropic satellite distribution and we consider 
the implications of our results, critically assessing the respective roles played by the 
orientation of the visible galaxy within its dark matter halo; anisotropic accretion 
of satellites from the larger scale environment; and the biased nature of satellites as 
tracers of the underlying dark matter subhalo population. 

Key words: galaxies: haloes — dark matter — galaxies: structure — galaxies: dwarf 
— galaxies: clusters: general — galaxies: formation 



1 INTRODUCTION 

The spatial distribution of satellites around isolated galaxies 
can provide important insights into the mass distribution in 
and around these galaxies. If dynamical effects could be ne- 
glected and if we could assume that satellite galaxies inhabit 
an unbiased set of dark matter subhaloes, then we would ex- 
pect satellites to cluster preferentially along the major axis 
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of the host dark matter hal o of the primary g alaxy, as in 
galaxy cluster mass haloes (|Knebe et alj |2004 ). We could 
therefore determine the orientation of the parent galaxy 
within its dark matter halo by using the spatial distribu- 
tion of its satellite galaxies. 

However, there is good reason to believe that dynamical 
effects will play an important role in dete rmining the spa - 
tial distribution of satellites. For example, lHartwickl (|2000h 
has argued that the imprint of anisotropic infall of satel- 
lites along filaments is evident in the orbits of recently ac- 
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creted systems, and IPefiarrubia et all (|2002h have shown 
that the inclination of a satellite's orbit can determine its 
dynamical response to the disc or spheroid of its parent 
galaxy. There is also good reason to suspect that the re- 
lationship between satellites and the underlying dark mat- 
ter subhalo population may no t be straightforwa rd, as has 
been argued by, f or example, iGao et all l|2004 ) (but see 
IConrov et al] |2006). Indeed, current galaxy formation mod- 
els indicate that satellite galaxies represent a biased subset 
of subhaloes whose spatial distribution is very likely to be 
aniso tropic within their parent galaxy's dark m atter halo 
(e.g. IZentner et ail 120051 ; iLibeskind et al.l l2005h . although 
current high resolution models suggest that the orientation 
of the sat ellite system with resp ect to the dark matter halo 
is robust (|Libeskind et alj|2007h . 

These considerations suggest that interpretation of the 
spatial distribution of satellites may be more complex than 
we might naively expect. Yet, despite these complexities, 
satellite galaxies represent a powerful observational probe 
into the mass distribution around galaxies. Therefore, it is 
essential to determine robustly the spatial distribution of 
satellite galaxies, and to establish whether or not they show 
a preferential alignment relative to their parent galaxies. 
This is the aim of the current study. 

Locally we see strong evidence for the preferential 
alignment or anisotropic distribution of satellites relative 
to their primary galaxies. Both the Milky Way and M31 
have satellite populations that lie in great planes that 
ar e highly inclined t o their dis c s. Th is ha s been noted 
by iLvnden-Bell (|l976l ). lHartwickl (|2000T) and iKroupa et all 
il2005l) for the Milky Way; by iKoch fc Grebell l|2006h and 
iMcConnachie fc Irwinl [ 2006 ) for subsamples of satellites of 
M31; and bv lMetz et alj (|2007l ), who performed an analy- 
sis of the statistical significance of the planar distribution 
around both galaxies. The Milky Way and M31 are the only 
galaxies for which the sample of satellites is large enough 
that their spatial distribution is analysed directly (without 
requiring stacking to obtain a statistical sample). They are 
also the only systems for which we are certain that the satel- 
lites are physically associated with their primaries, and for 
which the three-dimensional positions of the satellites with 
respect to their primaries are known. In addition, proper 
motion measurements exist for a number of Milky Way satel- 
lites, which confirm that these systems are on polar orbits 
<|Palma et al.ll2002h . 

Analysing the spatial distribution of satellites in ex- 
ternal systems is generally more challenging because no 
more than one or two satellites are detected per primary 
galaxy, and because the three-dimensional location of the 
satellite with respect to its primary is uncertain. This mea- 
surement requires that primaries and their satellites be 
stacked to obtain statistical samples from which the pro- 
jected angular distribution of satellites about a "typical" 
pri mary is det e rmine d. Such an analysis was first performed 
by iHohnbergl (|l969l ) for 58 isolated spiral and lenticular, 
or late-type, galaxies and their 218 optical companions, of 
which 75 were expected to be physically associated. He 
found that satellites at projected radii less than 50 kpc 
were more often found near the poles (minor axis) of the 
primary's disc (Figure [1} • This preferential polar distribu- 
tion, referred to as the "Holmberg Effect" , was subsequently 



confirmed by IZaritskv et"al] (|l997bh (hereafter IZSFWl) at 
larger projected radii of 300 - 500 kpc using a sample of 115 
spectroscopically-confirmed satellite galaxies of 69 isolated 
late- types. 

The advent of large galaxy redshift surveys such as the 



2dF Galaxy Redshift Survey (2dFGRS; IColless et a! 



and the Sloan Digital Sky Survey (SDSS; lYork et al 



2001) 



2000) 



has enabled this issue to be revisited. These surveys provide 
an abundance of galaxies with spectroscopic redshifts, 
and several recent studies have sought to use numbers 
of galaxie s far i n exc ess of t hose that were available to 
IHolmberd (|l969l l and IZSFWl to address the question of 
spatial anisotro p y wit h significan tly im proved statistics. 
ISales fc Lambasl (|2004T ) (hereafter |SL04|) used over 2000 
satellites of almost 1500 primary galaxies in the 2dFGRS, 
and found that satellites around blue primaries tended 
to follow an isotropic distribution, whereas the locations 
of low- velocity (|A-u| < 160 km s _1 ) satellites around red 
primaries tended to align with the major-axis. This finding 
contr asts sharpl y with the p olar dist r ibution infe rr ed by 
both IHolmberd dl969l) and IZSFTAfl . lYang et all (120061) 
( hereafter I Y06h . I Azzaro et all (|2007l ) (he reafter lAPPZl) and 
lAgustsson fc Brainerdl i|2007T ) (hereafter IAB07I) d etected a 
similar anisotropy of satellites in the SDSS; 1y06| used over 
16000 groups of galaxi es selec ted to lie within the same 
dark matter halo, while lAPPZl and fABOTl studied satellites 
of isolated galaxies; intriguingly, the latter study found 
that the isotropic distribution around blue galaxies was 
composed of a major-axis alignment for satellites at small 
projected radii and a m inor-axis alignm ent for sat ellites 
at large projected radii. iBrainerdl (|2005h (hereafter |B05| ) 
also selected satellites of isolated galaxies in the SDSS 
(her largest sample contained approximately 3000 satellites 
around 2000 primaries) and found that they exhibited a 
major-axis distribu t ion. Similar results were obtain ed by 
iFaltenbacher et all i|2007l ). and IWang et al l (l2008r ), who 
examined me mbers of tens of tho usands of groups in the 
SDSS. Finally. IX^aro et all (|2006h (hereafter [AZPKf) found 
no evidence for anisotropy based on a smaller sample of 193 
satellites of 144 isolated late-type galaxies in the SDSS. 

These results, derived from both 2dFGRS and SDSS, 
would seem to suggest that the preferential alignment of 
satellites about t he po les of t heir primaries noted by both 
IHolmberd l|l969h and IZSFWl was a consequence of small 
nu mber st atistics. However, the spatial anisotropy reported 
bv lZSFWl was detected with a statistical confidence greater 
than 99%, suggesting that small number statistics were un- 
likely to be a problem. Furthermore, as we have noted al- 
ready, there is strong evidence for the preferential alignment 
of satellites about the pole of the Milky Way and evidence 
for subsamples of satellites lying in inclined discs around 
M31, systems for which we can be certain that the satel- 
lites are physically associated with their primaries and for 
which we know their spatial distribution. Both of these ob- 
servations raise the spectre that it is systematic rather than 
random errors that more strongly affect detection of spatial 
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Figure 1. (Left) According to the Holmberg Effect, satellites of 
disc galaxies tend to lie near the minor axis of the disc. This cor- 
responds to an angle between the major axis of the galaxy and 
the satellite (the "disc angle", 6) of greater than 45°. (Right) Se- 
lection cylinders around a potential primary galaxy, which lies at 
the centre of the cylinders. The outer isolation cylinder is marked 
with solid lines and has radius -Router and length 2 Pouter . The 
inner isolation cylinder is marked with dashed lines and has ra- 
dius -Ri nner and length 2i)j nner . The satellites are drawn from the 
shaded cylinder with radius i? S at and length 2 v Ba ± ■ 



anisotropies in the distribution of satellites around external 
galaxies. This leads us to the issue of sample selection. 

In the case of both the Milky Way and M31, we have 
an abundance of detailed information that allows us to state 
with confidence which galaxies can be considered satellites 
belonging to these hosts. In the case of external galaxies, we 
do not have such detailed information and so we must em- 
ploy selection criteria that allow us to identify which faint 
galaxy neighbours in projection are likely to be satellites of 
the primary. Differences arising from how satellite galaxies 
are selected will affect measurements of the angular distri- 
bution of satellites, and will therefore influence how these 
data are interpreted. 

To illustrate this point, we know that the member 
galaxies of groups and clusters tend to cluster about the 
major axis of the Brightest Gr oup or Cluster G alaxy (BGG 
and BCG respect i vely; e.g. iBinggelil 1 19821 ; IWestl 1 19891 ; 
iMandelbaum et al.l 120061 ). If we identify these group or 
cluster members as satellites of their primary galaxy, the 
BGG or BCG, then we would interpret this measurement 
as showing that group or cluster members preferentially 
align with their primary's major-axis. However, we do not 
expect groups or clusters to be dynamically relaxed and so 
satellites in these systems might not trace the potential of 
an equilibrium mass distribution. Therefore, the study of 
satellites around isolated galaxies must be highly successful 
at excluding such groups and clusters. The ability of 
selection criteria to identify the proper type of system and 
to suppress contamination is essential for the result to be 
considered robust and physical. 

Mock galaxy catalogues constructed from cosmological 
TV-body simulations provide a powerful method for assessing 
the reliability of selection criteria. Dark matter haloes are 
populated with mock galaxies following either a statistical 
app roach based on "Halo Occu pation Distributions" (HOD; 
e.g. iBerlind fc Weinberg! i 2002), or a physically motivated 
appr oach based on sem i-analytical galaxy formation models 
(e.g. ICole et ai]|2000l ). Both approaches are parameterised 
and so statistical properties of the mock galaxy population, 
such as the luminosity function and the 2-point correlation 



function, are fine-tuned to recover the observed properties 
of real galaxies in the 2dFGRS and the SDSS. This provides 
an ideal testbed for selection criteria. We adopt the "Condi- 
tional Luminosity Function" (CLF) formalism of I Yang et al.l 
(2003), a statistical approach based on HODs, to develop our 
suite of mock catalogues. The CLF formalism allows us to 
assign to each dark matter halo in a cosmological iV-body 
simulation a probability of hosting N galaxies with a total 
luminosity L, and a distribution of luminosities drawn from 
a Schechter function whose parameters depend on the halo 
mass M. Further details are presented in § 12.31 

We note that previous studies have used mock galaxy 
catalogues derived from t he CLF formalism to inves tigate 
the radial distribution (Ivan den Bosch et all I2005T ) and 
kinematics (Ivan den Bosch et al.ll2004l) of satellite galaxies. 
However, we use our mock catalogues to establish optimal 
selection criteria that preferentially pick out isolated 
systems of primary galaxies and their satellites. We explore 
the impact of different selection criteria on the nature of 
satellite samples (e.g. group or cluster members?) and to 
establish robust criteria that minimise the influence of 
interlopers and a primary's larger-scale environment. This 
allows us to assess selection criteria adopted in previous 
studies and to quantify the ability of these criteria to 
identify isolated systems. 

By taking such care in establishing the criteria by which 
our isolated systems are selected, we are able to address a 
number of important outstanding physical questions using 
a robust sample of galaxies taken from the SDSS. In par- 
ticular, we revisit the question of whether or not the an- 
gular distribution of satellite galaxies about their primary 
is anisotropic, and whether or not anisotropy is linked to 
the morphological type of the primary. Such a dependence 
would present an exciting possibility to study the connec- 
tion between subhaloes, haloes, and galaxy morphology. In 
§[5] we discuss the physical significance of our results in this 
context. 

We note also natural overlaps with studies that seek 
to measure the flattening of dark matter haloes. There 
is good evidence to suggest that the haloes of early-type 
galaxies are flattened along the minor axis of the light - 
based on studies of X- ray isophotes (e.glBuote et al]|2002l ) 
and weak lensing (e.g iHoekstra et al.l 120041 ') - whereas the 
results for disc galaxies are inconclusive. For example, the 
tidal stream of the Sagittarius Dwarf has be en used to 
infer that the Milky Way's halo is spherical (|lbata et alj 
l200ll ; iFellhauer et"atl I2O O6), flattened in the same sense 
as the disc ( Martmez-Delgado et alj |2004| ; I Johnston et al 



2005), or flattened in the opposite sense as the disc l|Helmi 



2004), depending on which part of the stream is modelled 



( Law et al.ll2005l ). We discuss the implications of our results 
for halo flattening in § [5] 

The structure of this paper is as follows: in §[2] we dis- 
cuss the sample selection criteria, and present our sample of 
satellites of isolated galaxies. In § [3] we present our results 
for the anisotropic distribution of satellite galaxies with re- 
spect to their primary galaxy, with respect to large scale 
filaments, and the relative alignment of isolated galaxies 
with large scale filaments. In § [4] we compare our results 
to previous studies and analyse how different selection cri- 
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teria provide constraints on the origin of satellite anisotropy. 
Finally, we discuss the interpretation of the results in terms 
of the shapes of dark matter haloes, anisotropic infall, and 
the formation of galactic discs in § [5] 



for the r < 18 galaxies that we consider. Angular diame- 
ter distances and distance moduli are calculated assuming 
Om = 0.3 and Q,a = 0.7 and are quoted in /i-independent 
units, where Ho = h 100 kms -1 kpc -1 . 



2 SAMPLE SELECTION 

The most critical part of the analysis is the process of se- 
lecting satellite galaxies of isolated primaries. We must si- 
multaneously: 

(i) Minimize the number of primaries that are not iso- 
lated, i.e. which do not dominate the dynamics of their en- 
vironment. 

(ii) Minimize the number of "interlopers", or satellite- 
primary pairs that do not represent physical satellites of 
the primary galaxy. 

(iii) Maximize the sample size, subject to the above con- 
straints. 

We use mock catalogues generated from cosmological 
simulations to refine our criteria to fulfil the above require- 
ments and to critically examine the selection criteria that 
have been used in previous studies. 



2.1 Observational Data 

2.1.1 Sloan Digital Sky Survey 

Our sample is drawn fr om the Sloan Digital Sky Survey Data 
Release 6 (SDSS DR6) (jAdelman- McCarthy et ai1l2008l ) . All 
primary survey objects classified as galaxies in the imaging 
data that satisfy the spectrosc opic targeting algor ithm of 



Strauss et ail 120021) or the 



Eisenstein et al] l200lh 



either the Main galaxy sample 
Luminous Red Galaxy sample 
considered. Only galaxies with valid spectroscopic redshifts 
(sciencePrimary = 1 and zConf > 0.85) that are also spec- 
troscopically classified as galaxies are considered as primary 
or satellite galaxies; however, galaxies without such spectra 
or which are spectroscopically classified as stars (which of- 
ten indicates that a foreground star appears near the centre 
of a galaxy in projection) are also used when evaluating the 
isolation criteria. We have excluded all objects with unreal- 
istic colours (differences between successive bands of at least 
5 magnitudes), as this always indicate a spurious detection 
at the magnitudes we consider. Petrosian magnitudes are 
us ed throughout , and are dereddened using the corrections 
in lSchleeel et al.1 (ll99Sf> and k-corrected using KCORRECT 
v4_l_4 (|Blanton et aljkoOol ) to the 01 r bancfl to minimise 
the effect of errors in the k-correction. The nominal sur- 
vey limit for spectroscopic targets in the SDSS Main galaxy 
sample is r < 17.77; however, the actual limit varies across 
the sky. Therefore, we conservatively treat it as only being 
complete to r < 17.5, although we make use of galaxies as 
faint as r — 17.77 when they are available. When we re- 
quire photometric redshifts, we use the Dl neural network 
photometric redshif ts available in the SDSS DR6 catalogue 
|Ovaizu et al]|2008j ). which are the most accurate available 
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0.1; for simplicity we use the 



2.1.2 Galaxy Classification 

It is important to separate spheroid-dominated early-type 
galaxies from disc-dominated late-type galaxies. The ma- 
jor axes of spheroidal galaxies are determined by their 
anisotropic velocity dispersions while those of disc galaxies 
are determined by their angular momentum; therefore, their 
orientation with respect to the dark matter halo may be dif- 
ferent. Furthermore, the dynamical effects of discs versus 
spheroids on satellite orbits may be different. Galactic ori- 
entations may also depend on their history, which is probed 
by their stellar populations. Indeed, many previous studies 
have found that the satellites around red and blue galaxies 
are distributed differently. 

We adopt the galaxy classification scheme of 
I Bailin fc Harris! 42008b), which is based on the global 
concentration of the light profile and the location of the 
galaxy on the colour-magnitude diagram. This method 
explicitly accounts for inclination effects, and has been vali- 
dated using high-q uality imaging dat a from the Millennium 
Galaxy Catalogue l)Liske et al] 20051 ). Galaxies congregate 
in three distinct regions of parameter space: Early-type 
galaxies are red, highly concentrated, and ellipsoidal; 
Intermediate-type galaxies are red, have intermediate 
concentrations, and contain discs; and Late-type galaxies 
are blue, have low concentrations, and are disc-dominated. 
Determining whether the satellite systems of Intermediate- 
type galaxies bear closer relation to those of Early- vs. 
Late-type galaxies will provide useful insights into their 
nature. 

An analysis of the results when other classification 
schemes are adopted is given in Appendix [C] We find that 
our qualitative results are unchanged for any reasonable 
method of splitting the sample, although the numerical mag- 
nitude of the effect can vary by ~ la depending on the clas- 
sification method. 



2.2 Definition of Selection Criteria 

The format of o ur selection criteria is based on 
iNorberg et al] l|2008l ) and is similar to that used in previ- 
ous studies; however, the details differ in several important 
ways. To be considered isolated, primaries must not have 
any comparably bright neighbours within a large surround- 
ing region, and must be much brighter than all potential 
satellites in the immediate vicinity. We define three cylin- 
ders around each potential primary (see Figure [IJ: 

• Outer isolation cylinder: All galaxies within a projected 
separation of AR < -Router and a velocity difference of 
\Av\ < Woutcr must be at least m outcr magnitudes fainter 
in r. 

• Inner isolation cylinder: All galaxies within a projected 
separation of A_R < Ri aaet and a velocity difference of 
\Av\ < Winner must be at least mi nner magnitudes fainter 

in r (m innC r > m utcr). 
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• Satellite cylinder: Satellites are galaxies within a pro- 
jected separation of AR < _R sat and a velocity difference of 
|Ai>| < Usat. Satellites must be at least m sa t magnitudes 
fainter in r. For our criteria, we enforce i? sat = -Rinncr, 

^T-sat — ^^inner , and V sa t — 2 dinner • 

The adopted values of the parameters are given in Table Q] 
To ensure that satellites are not associated with more than 
one primary (a situation we refer to as a "multi-homed" 
satellite) and that there are no near neighbours too lumi- 
nous to be satellites, it is important that -Rinner < ^Router, 

dinner = ^outcr; Rsat ^ dinner? ^sat ^ ^^inncr, and 77l sa t = 

fMinner- These sanity checks are not fulfilled by many of the 
criteria that have been used in previous studies. 
We also apply the following additional criteria: 

(i) All primaries must be at least Warmer magnitudes 
brighter than r = 17.5 to ensure that all potential bright 
neighbours are brighter than the local survey limit. 

(ii) The projected distance from the primary to the near- 
est edge of the photometric survey footprint must be at least 
-Router to ensure that any potential bright neighbours have 
been observed photometrically. The projected distance from 
the primary to the nearest edge of the spectroscopic survey 
footprint must be at least f? sa t to ensure that all potential 
satellites are equally likely to have been observed spectro- 
scopically and therefore that the survey edge does not im- 
pose an angular bias in the selected satellite population. 

(iii) Because of spectral incompleteness, ~ 10% of galax- 
ies that fulfil the requirements of the SDSS Main galaxy 
targeting criteria do not have observed redshifts. There- 
fore, there are potential primaries that would not be con- 
sidered isolated if it turned out that their non-spectroscopic 
neighbours are at the same redshift. This issue is partic- 
ularly important because the limited number of fibres per 
tile causes the fractional completeness to be lower in re- 
gions of high galactic density. There are a number of ways 
of treating such galaxies (hereafter referred to as "viola- 
tors"): one could assume that most do not lie at the same 
redshi f t as t he primary and simply ign ore their presence 
(|SL04l : IboI lAZPKt lAPPZt lYOrj ; IAB07I ). one could estab- 
lish a threshold such that if there are more than a number 
iVvioi then the chances that at least one is at the same red- 
shift as t he primary is high and therefore eliminate th ose 
primaries l|Norberg et al.ll2008l: [Herbert -Fort et al.ll2008l ). or 
one could eliminate all such primaries on the grounds that 
there is a chance that t hey are not isolated (equivalent to 
setting iV v i i = 0: IZSFWh . The existence of photometric red- 
shifts in the SDSS DR6 catalogue allows us to use a more 
sophisticated method of determining whether the violators 
are likely to be at the same redshift as the primary. We first 
query the NASA/IPAC Extragalactic Database (NEE0) for 
literature spectroscopic redshifts of the violators of all pri- 
maries that would otherwise be included in our sample. If 
no spectroscopic redshift is available then we consider the 
photometric redshift z v ioi, photo; if it is within 2<7 v i i, photo of 
the spectroscopic redshift of the primary, where o- vio i >p hoto 
is the estimated error on the photometric redshift in the 
catalogue, then we eliminate the primary. We ignore the 
presence of violators that do not exist in the photometric 

3 http:/ /nedwww. ipac.caltech.edu/ 



redshift database, as these are mostly galaxies that are not 
detected in one or more bands; such galaxies are highly un- 
likely to be truly bright physical neighbours of the relatively 
luminous nearby galaxies that constitute our sample of pri- 
mary galaxies. We have confirmed that excluding primaries 
with such neighbours does not alter our conclusions. 

(iv) The total luminosity of the satellites must not be 
more than / sa t times the luminosity of the primary, and sys- 
tems with more than A sa tmax satellites are discarded. This 
ensures that the primary galaxy dominates the satellite sys- 
tem. 

(v) Postage stamp images of each potential primary were 
examined by eye. 11 objects were removed, 4 of which suf- 
fered from catastrophically bad background subtraction due 
to nearby bright objects and 7 of which are major mergers 
in progress. 

Our choices for the relevant parameters are motivated 
by analysis using mock catalogues as described in § 12.3.41 
Where applicable, the values previous authors have used for 
the selection parameters are given in Table [T] For those cri- 
teria that use only one isolation cylinder, we characterize it 
as an "Inner" cylinder. We summarize the selection criteria 
use d by eac h previous study below: 

ZSFW used somewhat thinner and shorter cylinders 
than those used in this work, but the criteria fulfil the above 
sanity checks. They did not have redshifts and magnitudes 
for every galaxy in the field, but eliminated all primaries 
that appeared by eye to have potentially criteria-violating 
neighbours (i.e. A v i i = 0). There was no formal edge of 
the surveyed area, and primaries were chosen to be at least 
2.5 mag brighter than the POSS magnitude limit, so all po- 
tential bright neighbours were considered. No cut on the 
number or luminosity fraction of satellites was imposed, and 
only morphological late- types were included as primaries. 

ISL04I used only one isolation cylinder, which did not sat- 
isfy m sa t < m ou tcr- Therefore, relatively bright galaxies are 
allowed to be in the satellite region. Because of this defini- 
tion, satellites may be multi-homed. Only primary galaxies 
with absolute magnitude Aib, — 51og/i < — 18 were used. 
They did not impose any constraints on proximity to the 
survey magnitude limit, proximity to the survey edge, A v i i, 
or / S at, b ut used N s 

atmax — 

lB05l tested three selection c riteria ; her Sample 1 (SI) 
used isolation criteria based on |SL04| (though adopting a 
different value of the Hubble constant); her Sample 2 (S2) 
used one very wide isolation cylinder, which permits bright 
galaxies to be in the satellite region and satellites to be 
multi-ho med; and her Sample 3 (S3) used isolation criteria 
based on IZSFWl (though adopting a different value of the 
Hubble constant). There was no cut on proximity to the 
survey edge, A v i i, or A sa tmax in any of the samples. There 
was a cut o n / sa t of 1.0. 

lAZPKl used Sample 2 of iPrada et all (|2003l ). but 
restricted the primaries to be morphological late-types, 
adopted a redshift limit of cz < 11000 km s _1 , and only ex- 
amined prima ries w ith —20.5 < Mb < —19.5. These criteria 
are similar to lSL04l . except that because m sa t = m ou tcr, the 
satellite cylinder is not permitted to contain bright galax- 
ies and satellites can only belong to a single primary. Al- 
though they did not formally adopt a cut on proximity 
to the survey edge, they searched for bright neighbours in 
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Table 1. Parameters of Selection Criteria 



Parameter 


This Work 


ZSFW 


SL04 


SI 


B05 
S2 


S3 


AZPK 


APPZ 
SI S2 


AB07 


Outer isolation cylinder: 






















Pouter {h~ 1 kpc) 


1000 


750 








700 










V outer (km S _1 ) 


1500 


1000 








1000 












0.7 


0.7 








0.7 










Inner isolation cylinder: 
























500 


375 


700 


490 


2000 


350 


500 


i on 


OUU 


511 


u; nilpl . (km s ) 


1500 


1000 


1000 


1000 


1000 


1000 


1000 


1000 


1000 


1000 


dinner 


2.0 


2.2 


1.0 


1.0 


0.7 


2.2 


2.0 


1.0 


2.0 


1.0 


Satellite cylinder: 






















ftat {h- 1 kpc) 


500 


375 


500 


350 


350 


350 


350 


350 


350 


365 


tJ Ba t (km S _1 ) 


750 


500 


500 


500 


1000 


500 


500 


500 


500 


500 


m 8a t 


2.0 


2.2 


2.0 


2.0 


1.5 


2.2 


2.0 


2.0 


2.0 


2.0 




NED+photo-z e 


e 


















/sat 


0.2 






1.0 


1.0 


1.0 








1.0 


-^Vsatmax 


4 




4 
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Sanity checks: 






















Forbids multi-homed satellites" 


Yes 


Yes 


No 


No 


No 


Yes 


Yes 


No 


Yes 


No 


Forbids nearby non-satellites 6 


Yes 


Yes 


No 


No 


No 


Yes 


Yes 


No 


Yes 


No 


Avoids survey edge c 


Yes 


N/A 


No 


No 


No 


No 


Yes e 


No 


No 


Yes 


Avoids survey magnitude limit d 


Yes 


Yes 


No 


No 


No 


No 


Yes 


No 


Yes 


No 



"The criteria do not permit satellites to belong to more than one primary galaxy. 

^The criteria do not permit there to be bright (non-satellite) galaxies within the satellite cylinder. 

c The criteria do not permit primaries so near the edge of the survey that potential bright neighbours would lie outside the 
survey region. 

d The criteria do not permit primaries faint enough that potential bright neighbours fall below the local survey magnitude 
limit. 

e See text for further clarification. 



Ide Vaucouleurs et all (jl99ll ) (|RC3T l. which fulfils the same 
purpose. The narrow range of absolute magnitudes and red- 
shift limit serve the same purpose as our cut on galaxies near 
the survey magnitude limit. They did not make any cut on 

max - or f s at - 

lAPPZl tested two set s of c riteria. Their Sample 1 (SI) 
uses the same criteria as |B05| SI but without the cut on 
/ sat . Their Sample 2 (S2) is identical to lAZPKl except that 
a wider range of primary luminosities, —23 < M r < —21 is 
used, the redshift limit is extended to cz < 30000 km s" 1 , 
and no check for bright neighbours outside of the survey 
boundary is performed. 

IAB07I used criteria very similar to |B05| SI, but they 
adopted a slightly different value for the Hubble constant, re- 
stricted A'satmax to 9, and ensured that their primary galax- 
ies were not near the edge of the spectroscopic survey. 

iHolmberd (|l969h did not have redshifts for any of his 
galaxies, making it difficult to directly compare our s elec- 
ti on. It is also difficu l t to c ompare our selection with lYOfj 
or iFaltenbacher et alj (|2007l ). who used an iterated percola- 
tion algorithm rather than isolation criteria. Their selections 
were tuned using mock catalogues to minimize the num- 
ber of interlopers; however, they were not designed to find 
isolated galaxies and most of their systems should be con- 
sidered groups or clusters rather than satellite systems of 
isolated galaxies. 

Our criteria are more restrictive than other criteria that 
have been used to select isolated galaxies. The advantage of 
using such a large sample as SDSS DR6 is less the abil- 
ity to boost the statistics than the ability to be extremely 



conservative with our selection criteria and still retain an 
acceptable number of galaxies. Given that the disagreement 
between previous results is more likely due to systematic er- 
rors than statistics, such a rigorous treatment is essential to 
disentangling the nature of the disagreement and determin- 
ing the true distribution. 

2.3 Mock Catalogues 

2.3.1 The Conditional Luminosity Function <&(L\M) 

We construct our moc k galaxy cata l ogues following t he pre - 
scription presented in lYang et alj (120041) (hereafter | Y04I) . 
which is based on earlier s tudies bv lYang et alj (|2003l ) and 
Ivan den Bosch et al.1 l|2003h . This approach allows us to as- 
sign to each dark matter halo of mass M a probability of 
hosting a population of N galaxies of total luminosity L. 
This probability is governed by the "Conditional Luminosity 
Function" , $(L|M) , which is parameterised by the Schechter 
function, 

$(L\M)dL = S- fJ±) exp{-L/L*)dL. (1) 

The normalisation $*, characteristic luminosity L* , and 
faint-end slope a are all functions of halo mass M; ap- 
propriate expressions for these quantities are in AppendixlAl 

Given $(L|M), we compute various "observable" prop- 
erties of the galaxy population associated with an average 
dark matter halo of mass M: 
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Table 2. Cosmological ./V-body Simulations. Lb ox refers to the 
length of the simulation box; rrt par t is the particle mass, which 
depends on both Lbox an d the number of particles iVpart; e is 
the gravitational force softening; M m ; n is a (conservative) esti- 
mate of minimum halo mass that can be reliably resolved, based 
on convergence of the mass function; and L m j n is the minimum 
luminosity. 



-^box 


mpart 


€ 




-^min 


[A _1 Mpc] 




[h- 1 kpc] 




[h- 2 L Q ] 


35 


0.21 x 10 9 


2.7 


10 10 


1.1 x 10 8 


50 


0.62 x 10 9 


3.9 


3 x 10 10 


1.1 x 10 8 


70 


1.7 x 10 9 


5.5 


9 x 10 10 


1.1 x 10 9 


100 


4.96 x 10 9 


7.8 


25 x 10 10 


1.1 x 10 9 



• the mean number of galaxies (N) (M) (see equa- 
tion [Bl]); 

• the luminosities of the central galaxy L ccn and satellite 
galaxies L sa t (see equation IB3|l : 

• the morphological type of each galaxy (i.e. early- versus 
late- type; see Appendix [B]} . 

All of these properties can be recovered using analytic dark 
matter halo mass functions, but a cosmological A-body 
simulation is required to assign phase space coordinates to 
mock galaxies. In what follows, we describe briefly the main 
steps involved in constructing the mock catalogues. A more 
detailed description is provided in Appendix [Bl 



2.3.2 Populating Dark Matter Haloes with Galaxies 

We perform a series of cosmological A^-body simulations 
following the formation of structure in a ACDM model 
from 2=50 to 2=0. Each simulation contains 256 3 particles. 
We adopt cosmological parameters of f2 m =0.3, f2A=0.7 and 
h—0.7, and the power spectrum of i nitial density pertur- 
batio ns (generated using CMBFAST; ISeliak fc Zaldarriagal 
1996) is normalised assuming a mass variance of ag=0.9. 
Details of the runs are presented in Table [3] We quote re- 
sults using the five 100 Mpc boxes (labelled "A" through 
"E" ) , but have verified using the smaller boxes that this res- 
olution is sufficient to reproduce the statistical properties of 
the mock catalogues. 

These simulations provide the dark matter host haloes 
that we populate with mock galaxies. Haloes are identified 
at 2=0 using the friends-of-friends (FOF) algorithm with 
a linking length of 6=0.2 times the mean interparticle 
separation. For each halo we compute its virial mass and 
radius, and we record also the coordinates of its most bound 
particle and a list of all particles that reside within its virial 
radius. This is (in principle) all the information we need to 
construct our mock galaxy catalogues. 

We note that $*, L* , and a axe functions of halo mass 
M, but they also depend on the cosmological parameters. 
lYang et al.l (|2003l ) presented a number of different CLFs 
for different cosmologies and different assumptions regarding 
the free parameters, and we adopt those used by I Y04l . who 
assumed a flat ACDM cosmology with fi m =0.3, Qa=0.7, 
h=0.7 and a normalisation <rg=0.9. Precise values for the 
CLF parameters are given in Appendix [X] 



These parameters were chosen to recover the observed 
luminosity functions and correlation lengths of galaxies in 
the 2dFGRS (both as a function of their luminosity and 
their type), but we find that the clustering and luminosity 
functions of galaxies in our mock catalogues are in very 
good agreement with the corresponding SDSS correlation 
and luminosity functions. 

It is important to estimate the minimum reliably re- 
solved halo mass M m j n ; below this threshold the number 
density of haloes tends to be suppressed relative to the num- 
ber density they would have in the limit of infinite numerical 
resolution. In this work we assume that the mass function 
is converged for haloes containing 50 particles or more (see 
discussion in Appendix[B}; this gives M m ; n =50 m par t in Ta- 
bled 

Knowing M m i n allows us to estimate .Lmin, the mini- 
mum luminosity that we can assign to a mock galaxy, us- 
ing the "conditional probability distribution" P(M\L). L m i n 
fixes the halo occupation number, the average number of 
galaxies per halo of mass M, {N} (M) (equation IB1|) ; as 
M m i n and therefore I/ m i n increases (decreases), (N) (M) in- 
creases (decreases). The mass resolution of our simulations 
means that we adopt L m i n = 1 .1 X 10 9 h~ 2 Lq in our mock 
catalogues A to E, upon which our analysis is based. 

We assume that the number of galaxies A^ in a dark mat- 
ter h alo is Poisson distributed about (N) (M) (|Yang et al.l 
2003), and that each galaxy is assigned a luminosity drawn 
from <f>(L|M). The central galaxy is defined to be the bright- 
est galaxy in the halo and has a luminosity L > L\, where 
L\ satisfies the condition that (N) (M) = 1 (see equa- 
tion lB3|) . The remaining N-l galaxies have luminosities that 
are drawn randomly from the luminosity function in the 
range L min < L < Li. 

The morphological type of each galaxy is obtained from 
/iate(i, M) (equation IB4|) . the probability that a galaxy of 
luminosity L hosted by a dark matter halo of mass M is late- 
type. Finally, the position and velocity of the central galaxy 
are associated with those of the most bound particle in the 
halo, while the positions and velocities of the remaining iV-1 
galaxies are obtained by randomly sampling the particles in 
the FOF group. 

A more detailed description of our approach to populat- 
ing dark matter haloes with galaxies is given in Appendix iBl 

2.3.3 Constructing a Mock Galaxy Redshift Survey 

At this point we depart from |Y04| . who wished to study 
a mock 2dFGRS and stacked simulation boxes to recover 
the survey's median redshift. Our needs are more modest 
- we wish to evaluate the reliability of our sample selection 
criteria. To transform our raw mock galaxy distribution into 
a mock galaxy catalogue, 

• We select a single simulation box, recalling that each 
box has periodic boundary conditions, and we replicate it 
3 times along each dimension, producing a stack of 3x3x3 
boxes. We then centre the stacked boxes on the median red- 
shift of the SDSS 2 mcd =0.11. 

• We place a virtual observer at 2=0 and we define an 
(a, S) coordinate frame with respect to the centre of the 
stacked box. 
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• We compute a redshift z for each galaxy as seen by the 
virtual observer from the recessional velocity cz = Hr + v ■ 
f/\r\, where the galaxy is at r with respect to the virtual 
observer, and v-f/\r\ is its line-of-sight velocity. We account 
for observational velocity uncertainties by adding a random 
velocity drawn from a Gaussian of width 30 km s _1 . We also 
compute the apparent magnitude according to its luminosity 
and distance, to which we add an RMS error of 0.02 mag, in 
accordance with the SDSS internal estimates of the redshift 
and photometric errors for galaxies. 

• We remove the redshifts of all galaxies that fall below 
the magnitude limit r = 17.77 of the SDSS Main galaxy 
sample and those with declinations less than +30° , for which 
the depth of the stacked simulation boxes is insufficient and 
which we define as the edge of the mock survey boundary 

There are a number of attributes of the observed survey 
that are not accurately reproduced by the mock survey. The 
true survey boundary is much more complicated than the 
boundary of our mock catalogue, and hence a much larger 
fraction of galaxies lie near an edge and may be near an 
unseen bright galaxy. Also, the true survey is not spectro- 
scopically complete, and is less complete in regions of higher 
density due to the lower availability of fibres. Therefore, our 
analysis using the mock catalogues would underestimate the 
importance of excluding primaries near the survey edge, ex- 
cluding primaries near the magnitude limit of the survey, 
and implementing a cut based on N v - 10 i. As these issues can- 
not be addressed well using the mock catalogues, we do not 
implement them in our analysis of the mock catalogue. Their 
effects on the sample are investigated empirically in § 13.11 



2.3.4 Tests of our Selection Criteria 

Using the mock catalogues, we quantify the degree to which 
current and previous selection criteria accurately identify 
physical satellites of isolated galaxies. We quote the results 
of the "Mock A" catalogue (see Table ID1[I , but the results 
are consistent with those from the other mock catalogues. 

Interlopers are identified in the mock catalogues as 
satellites that do not belong to the same halo as their se- 
lected primary galaxyQ We do not consider primary galaxies 
as isolated if: 

(i) they are not the central galaxy of their halo, or 

(ii) they are not sufficiently more massive than other 
galaxies within the halo. 



Point (i) is easily determined from the mock catalogues be- 
cause we know which galaxy is the central galaxy of each 
halo. Point (ii) is more difficult to determine because the 
CLF formalism assigns luminosities and not masses to in- 
dividual galaxies. However, we can estimate the degree to 
which the primary galaxy dominates its halo by calculating 
the fraction of the total halo luminosity that is contributed 
by the primary galaxy, 



4 Note that many previous works refer to all unwanted satellites 
as "interlopers", regardless of the reason for wanting to exclude 
them from the sample. We prefer to designate only satellites that 
are not physically associated with their selected primary as "in- 
terlopers" to distinguish them from satellites that are physically 
associated with unwanted primaries. 




Figure 2. (Top) Histogram of the number of primaries selected 
from the mock catalogue as a function of / pr i m , the fraction of 
the true halo luminosity that comes from the primary galaxy. The 
solid lines represent systems chosen using our adopted selection 
criteria, while the other line styles indicate othe r re prese ntative 
criteria (dashed, dotted and dot-dashed for |B05| S2, bos! S3 and 
lAPPZl S2). (Bottom) As above, but weighted per satellite in the 
sample. The vertical dotted lines denote the / pr i m below which 
primaries are considered non-dominant. 



/prim — Tpyi m /Lt- 



(2) 



In the top panel of Figure [2] we plot a histogram of the 
number of selected primarie s as a fu nctio n of /p rim for f our 
sample sets of criteria (ours. lB05l S2. IB05I S3 and lAPP^ S2) 
that span the range of observed behaviours. In the bottom 
panel, we weight each primary by its number of selected 
satellites to demonstrate its influence on the observed sam- 
ple. Because primaries that contribute less to the luminosity 
of their halo have more satellites, the tail to low fractions 
is exacerbated. The histograms are characterised by a sym- 
metric peak centered at / pr im = 0.9 that extends down to 
0.8 that we identify as truly isolated primaries, and a long 
tail to low values that we wish to eliminate. Based on ex- 
amination of these histograms, we consider a primary to be 
"non-dominant" if / pr i m < 0.8 (denoted by the vertical line 
in Figure [2j|. 

Table [3] lists the number of primaries and satellites se- 
lected from the mock catalogue that pass and fail the above 
interloper and isolation criteria. Row 1 indicates the number 
of primaries selected, row 2 indicates the fraction of those 
primaries that are not the central galaxy of their halo (cat- 
egory (i) above), row 3 indicates the fraction of primaries 



that do not dominate the dynamics of their halo (category 
I (ii) I above), and row 4 indicates the fraction of selected pri- 
maries that are truly isolated. Row 5 indicates the number 
of satellites selected, row 6 indicates the fraction of those 
satellites that are "multi-homed" , i.e. selected as satellites of 
more than one primary, row 7 indicates the fraction of satel- 
lites that do not belong to the same halo as their selected 
primary ( "interlopers" ) , row 8 indicates the fraction of satel- 
lites that are physically associated with non-dominant pri- 
maries ("/non-dom"), and row 9 indicates the fraction of 
selected satellites that are selected correctly, i.e. they are 
physically associated with isolated primaries. 

When evaluating criteria based on different bandpasses, 
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we use the following simple t ransfo rmations: for the AZPK 
criteria, which are based on IRC3I Bt magnitudes, we as- 
sume a constant B — 01 r = 0.6, ty pical of the late-type 
galaxies they studied; for the lAPPZl criteria, we assume a 
constant difference between the z — and z = 0.1 r-band k- 
corrections of 0.23 mag; and for the lSLoj criteria, which are 
based on UK Schmidt bj magnitudes, we assume a constant 
bj -°'V = 0.7, intermediate between the typical values for 
early and late type galaxies. 

The use of the mock catalogues to evaluate selection 
criteria for studies that did not use SDSS data is not as 
accurate as for those studies based on SDSS because the 
parameters of the mock catalogues, such as the typical pho- 
tometric and redshift error, are sp ecificall y tuned to mimic 
SDSS. T he ph otometric errors in ZSFW and the velocity 
errors m SL04 are significantly larger. 

All sets of criteria do an adequate job of selecting cen- 
tral galaxies as primaries and of minimizing the fraction of 
interlopers and multi-homed satellites, with each source of 
contamination contributing less than 5% to each sample- 
However, the fraction of the sample that lies around non- 
dominant primaries (/non-dom) extends from a low of 6% for 
our adopted criteria to almost 85% in the case of |B05| S2. 
We have examined by eye the fields surrounding a subset of 
the primary galaxies in our SDSS sample in order to confirm 
that this is an accurate measure of the degree of isolation of 
the sample and estimate that 7% of our primaries are mem- 
bers of groups, in very good agreement with the value that 
we derive from the mock catalogues. 

The greatest single predictor of the magnitude of 
/non-dom is whether non-satellites are permitted to lie within 
the satellite cylinder. In other words, criteria with m sat > 
JTiinner generally fail to select isolated primaries. Although 
we do not evaluate the effects of iVvioi using the mock cat- 
alogues, neglecting to account for spectroscopic incomplete- 
ness is the other factor that can result in galaxies larger 
than satellites lying within the satellite cylinder; we there- 
fore expect that this also has a significant effect on /non-dom- 
The impact of non-dominant primaries on the results will be 
discussed in detail in §|4] 

One might ask how sensitive these conclusions are to the 
CLF method used to assign luminous galaxies to dark mat- 
ter halos versus, for example, a semi-analytic model. Moti- 
vated by the halo luminosity funct ion found i n a pa rticularly 
discrepant semi-analytic model of lEke et al.l (|2004l . their fig- 
ure 5), we have tested a mock catalogue where we arbitrarily 
doubled the luminosity of the central galaxy. The resulting 
/non-dom changes by less than 0.1 for the vast majority of 
selection criteria tested. Given the insensitivity of / n0 n-dom 
to such relatively dramatic departures from our method of 
assigning luminous galaxies to dark matter halos, we feel 
confident that our conclusions regarding the selection crite- 
ria used in previous studies are robust. 

We have optimized the parameters we use for the se- 
lection criteria using the mock catalogues to maximise both 
the size of the sample (Row 5) and the fraction of the sample 
that passes all of the checks (Row 9). In particular, all pa- 



5 Although it is in principle possible to select multi- homed satel- 
lites using many of the sets of criteria, only for lB05l S2 does this 
situation ever occur in practice. 



rameters in Table [T] were varied in turn and the new value 
kept if it increased the sample size without a correspond- 
ingly large increase in either the incorrectly-selected fraction 
(the inverse of the value in Row 9), the interloper fraction, 
or the fraction of satellites around non-dominant primaries; 
or, conversely, if it decreased the incorrectly-selected frac- 
tion without a correspondingly large decrease in sample size. 
When in doubt, we erred on the side of more restrictive cri- 
teria. This process was repeated until convergence. 

2.4 The Sample 

The following quality cuts are imposed on the sample: 

(i) Satellites within 35/i _1 kpc of the primary are re- 
moved due to the known sky subtraction pro blem around 
bright sources (iMandelbaum et al. 1 120051 . 120061 ) and the pos- 
sibility of bright knots in the outer regions of a galaxy being 
mistakenly deblended as separate galaxies. 

(ii) As the interloper fraction rises at large projected 
radius, and the radius at which interlopers dominate in- 
creases with halo mass, we eliminate satellites of less lu- 
minous [M r — 5 log ft > —21.1, i.e. fainter than the me- 
dian) intermediate- and early-type primaries and of all 
late-type primaries which lie beyond a projected radius of 
345 kpc; this choice of parameters is justified below. 

The resulting sample of primary and satellite galaxies 
is given in Table [4] The full sample contains 866 satellites of 
722 primaries; 311/138/273 of the primaries are classified as 
ear ly/intermediate/late- type hosting 378/167/321 satellites. 

The following cuts are further imposed on systems used 
to measure the anisotropy around primary galaxies: 

(i) Primaries that do not have a measured position an- 
gle (PA) in the SDSS DR6 database are excluded; this cut 
excludes 11 primaries. 

(ii) Primaries with isophotal axis ratios b/a > 0.8 are ex- 
cluded. Galaxies with nearly circular isophotes have poorly- 
constrained PAs. In addition, any anisotropy that exists in 
three dimensions gets washed out in projection as the sys- 
tem is viewed close to its axis of symmetry. The numerical 
choice of b/a > 0.8 for this cutoff is motivated in § 13.11 and 
the effects of changing this value are discussed. 

Unless otherwise specified, this subsample is the sample re- 
ferred to for the remainder of the paper, and contains 440 
satellites of 372 primaries. In § 13.21 we analyze the distri- 
bution of satellites with respect to the large scale structure 
(LSS). For those purposes, the following further cuts are 
imposed instead: 

(i) Systems within a projected radius 3000 h~ kpc of the 
edge of the spectroscopic survey footprint are excluded in 
order to ensure that the survey boundary does not introduce 
a bias in the derived LSS axis. 

(ii) Primaries for which the LSS axis is undefined because 
there are no galaxies within the cylinder used to define the 
axis or for which the LSS axis ratio is greater than 0.9 are 
excluded. 

This subsample contains 572 satellites of 480 primaries. 

The distribution of luminosities, number of satellites 
per primary, radial separations and velocity differences are 
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Table 3. Results of Applying Selection Criteria to Mock Catalogues 



Parameter 


This Work 


ZSFW 


SL04 




B05 




AZPK 


APPZ 


AB07 










SI 


S2 


S3 




SI 


S2 




Selected primaries 


337 


135 


1516 


1828 


1778 


404 


114 


1595 


636 


589 


Non-central primary fraction 


0.000 


0.007 


0.001 


0.002 


0.002 


0.002 


0.000 


0.003 


0.000 


0.005 


Non-dominant primary fraction 


0.071 


0.059 


0.600 


0.640 


0.724 


0.087 


0.123 


0.656 


0.272 


0.626 


Isolated primary fraction 


0.929 


0.941 


0.400 


0.360 


0.276 


0.913 


0.877 


0.344 


0.728 


0.374 


Selected satellites 


388 


187 


2396 


3461 


3852 


563 


133 


3081 


980 


1112 


Multi-homed fraction 


0.000 


0.000 


0.000 


0.000 


0.001 


0.000 


0.000 


0.000 


0.000 


0.000 


Interloper fraction 


0.046 


0.048 


0.037 


0.036 


0.038 


0.025 


0.023 


0.042 


0.036 


0.049 


/non-dom 


0.062 


0.102 


0.691 


0.767 


0.844 


0.151 


0.218 


0.778 


0.409 


0.754 


Correctly-selected satellite fraction 


0.892 


0.861 


0.290 


0.216 


0.138 


0.824 


0.759 


0.203 


0.563 


0.224 



shown in Figure [3] In Figure [3^ we plot the absolute mag- 
nitude distributions of the primary and satellite galaxies, 
along with the magnitude difference. The median absolute 
magnitude is M r — 51og/i = —21.1. The typical primary 
has a luminosity similar to the Milky Way, while the typ- 
ical satellite is ~ 2.5 magnitudes fainter, slightly brighter 
than the LMC. Figure [3]d shows the number of satellites per 
primary. Most primaries are surrounded by only one satel- 
lite. Figure presents the distribution of radial separation 
between satellites and primary galaxies for samples split by 
the type and luminosity of the primary (the median absolute 
magnitude, —21.1 is used as the cutoff between the "bright" 
and "faint" subsamples). In this panel only, we include the 
distant satellites of faint primaries. The radial distribution 
of true satellites declines with radiu s, while the distrib ution 
of interlopers increases with radius (jChen et all [2006:1. The 
crossover between these regimes occurs at the edge of the 
dark matter halo, and therefore depends on galaxy mass. 
Figure |3j; shows no evidence for an increase in the number 
of outer satellites due to interlopers around bright early- 
and intermediate-type primaries. However, such an upturn 
is evident around faint primaries of both morphologies and 
all late- type primaries beyond AR > 345 h~ l kpc. There- 
fore, around these primaries we exclude satellites at pro- 
jected radii AR > 345 h' 1 kpc in all other panels and re- 
maining plots; however, we have confirmed that if we include 
these satellites, our results are qualitatively unchanged. Fig- 
ure [3}I presents the distribution of velocity differences be- 
tween satellites and their primaries. The velocities of se- 
lected satellites cluster strongly about the velocity of the pri- 
mary indicating that they are indeed physically associated. 
The velocity dispersion of satellites around bright primaries 
is higher than around faint primaries due to their la rger 
mass (see, e.g. IConrov et al.ll2007l : iNorberg et al.ll2"o"osl . and 
references therein). 

3 RESULTS 

3.1 Distribution About The Primary Galaxy 

We define the "disc angle" as the angle between the r-band 
isophotal position angle of the major axis of the primary and 
the PA of the great circle between the primary and satellite 
(see Figure [TJ; we fold this angle into the range 0-90°. If 
satellites are distributed isotropically around the primaries 



then the distribution of disc angles is uniform with a mean 
of 45°. 

In Figure [4] we plot the cumulative and differential disc 
angle distributions. The full sample shows a tendency to lie 
at small disc angles, i.e. for the satellites to lie near the major 
axis of their parent galaxy. The hypothesis that the angles 
are distributed randomly is ruled out at a greater than 2a 
level: according to the Kolmogorov-Smirnov (KS) test, the 
probability is 0.06 with a mean disc angle of 42?5 ± 1?2. 

The satellites surrounding different types of primaries 
show different angular distributions: there is a clear excess 
concentration of satellites along the major axis for early- 
types, hints of a major-axis excess around intermediate- 
types (the measured magnitude of the anisotropy is in fact 
larger than around early-types, but is detected at less than 
2a), and no detectable anisotropy around late- types. The 
mean disc angles are 41? 3 ± 0?9, 40? 5 ± 2? 5 and 45? 5 ± 2?1 
around early, intermediate and late-types respectively, with 
KS test probabilities of being drawn from a random distribu- 
tion of 0.03, 0.10, and 0.99. When the samples are combined 
these effects counteract each other, with the isotropic satel- 
lites of late-types diluting the major-axis alignment seen 
around the early-types. The KS test results, the mean and 
median disc angles, and the polar fraction (the fraction of 
satellites with disc angles larger than 45°) are listed in Ta- 
ble [5] The quoted uncertainties are determined using boot- 
strap resampling of the primaries and represent 68% confi- 
dence intervals. All statistics reinforce the same conclusion: 
satellites of spheroidal early-type galaxies tend to lie near 
the long axis of the spheroid, while satellites of disc galax- 
ies are isotropically distributed around blue late-types but 
show hints of lying near the disc around red intermediate- 
types. We have verified using the mock catalogues, which 
have intrinsically isotropically-distributed satellites, that we 
could not have measured this level of anisotropy around the 
early- type primaries if it were not physically present; see 
Appendix [D] 

Our measured alignment must be a lower limit on the 
intrinsic three-dimensional alignment for several reasons. 
Firstly, images of triaxial elliptical galaxies contain isopho- 
tal twists due to projection effects, which introduces scatter 
between the isophotal PA we use and the intrinsic 3D major 
axis. Secondly, galaxies seen closer to their symmetry axis 
have their anisotropic signal diluted. Finally, interlopers may 
dilute the signal. 

In Figure [5^, we plot the mean disc angle of satel- 
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Figure 3. (a) Absolute magnitude distributions of primary galaxies (thick lines) and of satellites (thin lines) (left panel) and the difference 
between the magnitude of each satellite and its primary (right panel). The red/dotted lines shows the distribution for early-type primaries, 
the green/dot-dashed lines shows the distribution for intermediate- type primaries, and blue/dashed lines refer to late- type primaries. Note 
that the satellites are also separated according to the type of their primary, not according to their own type, (b) Distribution of number 
of satellites per primary. Line colours/styles indicate the classification of the primary galaxy as in (a), (c) Distribution of projected radial 
separation of satellites from their primary. Left panels show satellites of early-type galaxies, middle panels show satellites of intermediate- 
type galaxies, and right panels show satellites of late-type galaxies. Top panels show satellites of primary galaxies more luminous than 
the median [M r — 5 log h < —21.1) while bottom panels show satellites of the less-luminous primary galaxies (M r — 5 log h > —21.1). The 
vertical lines mark the maximum radial extent for satellites. These satellites are not included in the remaining plots, (d) Distribution of 
velocity differences between satellites and their primary. Panels are as in (c). 



Table 5. Anisotropy of Satellite and Primary Distributions 



Parameter 



Full sample Early-type primaries Intermediate-type primaries Late-type primaries 



Disc angle: 

KS probability 

Mean disc angle [°] 

Median disc angle [°] 

Polar fraction 

Early-type satellite mean 

Intcrmcdiate-type satellite mean 

Late-type satellite mean 
LSS angle: 

KS probability 

Mean LSS angle [°] 
Disc/LSS angle: 

With satellites - KS probability 

- mean 

All isolated — KS probability 

- mean 



0.06 

42.5 ± 1.2 
40.9li;9 

0.46 ±0.02 

42.5111 
40.8±37 

42.8 ± 1.3 

0.56 

44.6 ± 1.2 

0.06 
41.8 ± 1.4 

0.0001 
42.6 ±0.7 



0.03 
41.3 ±0.9 
37.1±|;| 
0.44 ±0.03 



36.6 



+6.0 



6.7 
,4-4.0 
-4.5 

41.6 ±2.3 



42. r 



44.0 



0.76 
+ 1.5 
1.7 



0.4 
43.4 ±2.2 

0.01 
42.1 ± 1.2 



0.10 
40.5 ±2.5 
40.5^-1 

0.44 ±0.05 



49.5 
43.0 
38.9 



+6.6 
7.4 
+7.9 
7.6 
+2.6 
2.8 



0.49 



43.0 



+3.0 



0.05 



39.8+22 
0.0006 
41.3 ± 1.2 



0.99 

45.5 ±2.1 
43.6+i;? 

0.50 ±0.04 

47 7+11-7 
^'■'-10.6 

28.7±S;| 
46.2 ±2.1 

0.40 
46.4 ±2.2 

0.10 
41.7 ±2.3 
0.17 

43.6 ±0.9 



Satellite Galaxy Anisotropy 13 




15 30 45 60 75 90 15 30 45 60 75 90 
Disk angle [degrees] Disk angle [degrees] 

Figure 4. (Left) Cumulative distribution of angle between the major axis of the primary and the location of the satellite ("disc angle"). 
The thick black/solid, red/dotted, green/dot-dashed, and blue/dashed lines refer to the distribution of satellites around all primaries, 
early-type primaries, intermediate-type primaries, and late-type primaries respectively. The thin dotted line shows the distribution 
expected if satellites were isotropically distributed. (Right) Differential distribution of the disc angles. The error bars are determined by 
bootstrap resampling of all primary galaxies. 



lites as a function of projected separation from their pri- 
mary. Satellites are binned in 3 annuli spaced evenly in ra- 
dius from 35 ft -1 kpc to 500 kpc, and are separated into 
those more luminous and less luminous than the median 
M r — 5 log ft = —21.1. There is no significant trend with ra- 
dius. Our intermediate-separation bin covers approximately 
the same radii as the outermost bin of IAB07I . in which they 
detect minor-axis alignment around late-type primaries; our 
mean disc angle around late-types at these radii is also 
greater than 45°, but not at a statistically significant level. 
While we do not detect the major-axis alignment that they 
see at small radii, a large number of systems in their in- 
nermost radial bin fall within the 35 h~ l kpc region that 
we exclude to avoid contamination from Hn regions in the 
outskirts of the parent galaxy; if such regions are mistak- 
enly included as satellites, they will bias the result towards 
major- axis alignment. 

In Figure [5]d, we compare the anisotropy as a function 
of the luminosity of the primary, luminosity of the satellite, 
and of the difference in magnitude between the primary and 
satellite. The widths of the bins are chosen such that there 
are similar numbers of satellites in each bin; in all panels, 
the symbols are plotted at the mean luminosity or mean 
AM r of the galaxies in the bin. No clear trend is apparent 
as a function of primary luminosity. The anisotropy shows 
a weak dependence on the magnitude difference between 
the primary and satellite, although the sense of the trend 



differs for different primary types. Around early-type pri- 
maries, those satellites that are not much fainter than their 
primaries show a stronger major-axis alignment than the 
satellites that are much fainter, while around intermediate- 
type primaries the opposite trend holds. Around late-type 
primaries, those satellites that are bright relative to their 
primary show ~ 2a major-axis alignment, while the faintest 
50% of the satellites relative to their primary do not (or, if 
anything, show a polar distribution, although not at a statis- 
tically significant level). These trends could either reflect a 
dependence on the degree to which the primary dominates, 
i.e. it could truly depend on the magnitude difference, or it 
could reflect a dependence on the luminosity of the satellite. 
The right panel of Figure 0) reveals that the anisotropy 
around intermediate-type galaxies can equally well be ex- 
plained as being a function of satellite luminosity, while the 
angular distributions around early- and late-type primaries 
are not; therefore, for these galaxies, the trends seen in the 
middle panel truly reflect a dependence on the relative dom- 
inance of the primary. 

IKoch fc Grebell l|200rj ) found that only t he e arly- 
type satellites of M31 h ave a polar alignment. IyoI and 
iFaltenbacher et "aD (|2007l ) found that the red satellites of 
red primaries sho w stron ger major-axis alignment than the 
blue satellites, an d SL04 found that quiescent satellites show 
stronger anisotropy than star- forming satellites. To deter- 
mine if these signals are evident in our sample, we split the 
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Figure 5. (a) Mean disc angle of satellites as a function of pro- 
jected separation from their primary. Black plusscs, red diamonds, 
green trianges and blue squares refer to satellites of all primaries, 
early-type primaries, intermediate-type primaries, and late-type 
primaries respectively. Radial bins are 155 h~ 1 kpc wide and are 
plotted at the central bin radius with different symbols offset 
slightly for clarity. The top panel contains all satellites while the 
middle and bottom panels contain satellites of primaries with 
M r - 51og/i < -21.1 and M r - 51og/i > —21.1 respectively, (b) 
Mean disc angle of satellites as a function of the absolute mag- 
nitude of the primary (left), the magnitude difference between 
the primary and satellite (middle), and the absolute magnitude 
of the satellite (right). Bins are chosen to have approximately 
equal number of satellites per bin and are plotted at the mean 
magnitude. 

sample by the galaxy type of both the primary and satellite 
and plot the distributions in Figure [6] Mean disc angles as 
a function of satellite type are given in Rows 5 through 7 
of Table [5] Although we see small deviations for particu- 
lar subsamples (e.g. the early-type vs. late-type satellites of 
intermediate-type primaries, or the intermediate-type satel- 
lites of late- type primaries) , none are statistically significant 
due to the small number of early- and intermediate-type 
satellites. 

In § 12.41 we excluded primary galaxies with isophotal 
axis ratios b/a > 0.8. Our choice of cutoff is motivated by 
Figure [7J where we plot the mean disc angle around pri- 
maries with axis ratios less than or equal to the plotted ab- 
scissa. A cutoff at b/a = 0.8 provides enough systems that 
the results have converged, while avoiding the dilution in the 
signal seen at higher values of b/a. For a circular disc galaxy 
of intrinsic thickness 0.2, this corresponds to an inclination 
of 43°. 

As discussed in § 12.3.31 it is difficult to assess the effects 
of some of the selection parameters using the mock cata- 
logues. Therefore, we now empirically investigate the effects 
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Figure 6. Cumulative (top) and differential (bottom) distribu- 
tions of disc angle as a function of morphological type of both 
the primary and satellite. The three panels specify the type of 
the primary, while line colours and styles denote the type of the 
satellite. 
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Figure 7. Mean disc angle for satellites of primaries with isopho- 
tal axis ratios less than or equal to the plotted abscissa. One error 
bar is plotted at the location of each primary galaxy. Colors/line 
styles are as in Figure [4] 

of varying all of the selection parameters, including several 
parameters that have been neglected by many previous stud- 
ies. In Figure [8] we plot the mean disc angle and the sample 
size as we adjust N v - 10 \, / sa t, and whether or not the survey 
magnitude limit or survey edge are taken into account. The 
mean disc angles are shown as symbols with error bars while 
the early and late-type sample sizes are shown as histograms 
above and below (the intermediate-type sample sizes show 
identical trends). 

In panel (a), we compare our combined NED and pho- 
tometric redshift-based method of dealing with criteria vio- 
lators, a method that only uses the photometric redshift, 
and methods based on a cut at various values of N v i \. 
More restrictive values of JVvioi lead to smaller sample sizes, 
particularly for iV v i i < 4. The measured anisotropy is 
relatively constant as a function of N v - lo \. Our combined 
NED+photometric redshift method produces a sample size 
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equivalent to using iV v j i = 1, and therefore recovers a rea- 
sonable fraction of the sample available from ignoring the 
presence of violators, while conservatively excluding any sys- 
tems that could contaminate the sample. In panel (b), we 
plot the effects of / sat . The measured anisotropy rises as / sat 
is reduced to very low values; this is accompanied by a sharp 
decrease in sample size. As indicated by Figure [5p, this in- 
dicates that in many cases the anisotropy is related to the 
dominance of the primary galaxy, and samples selected with 
large values of / sat may be significantly contaminated. In 
panel (c), we show the effects of ignoring the survey magni- 
tude limit or of ignoring the survey edge. Both samples are 
marginally larger with no significant change in the measured 
anisotropy. 

Although the exact strength of the anisotropy can vary 
with some of these oft-neglected parameters, the qualitative 
result, that the satellites of early- and possibly intermediate- 
type galaxies show a major-axis distribution while the satel- 
lites of late-type galaxies are isotropically distributed, is not 
dependent on the value of any one of these parameters. Vary- 
ing the other parameters from Table [1] has very little effect 
on t he measured anisotropy, as anticipated by the results of 
|B05| . who found identical results in three samples with quite 
different values of these parameters. 



3.2 Satellite and Primary Distribution Relative to 
Large Scale Structure 

If satellites are accreted from filaments, then the most 
recently-accreted satellites will be aligned preferentially with 
the surrounding filamentary large scale structure (LSS). We 
test this expectation by determining the axis of the LSS sur- 
rounding each primary galaxy. To determine this axis, we use 
all spectroscopic galaxies with projected radii between 1000 
and 3000/i -1 kpc (thereby explicitly ensuring that there is 
no overlap between the galaxies used to determine the ori- 
entation of the LSS and those used to evaluate the isolated- 
ness of the primary or the satellites themselves) with veloc- 
ities that differ from that of the primary by no more than 
400 km s . The velocity dimension of this cylinder is sig- 
nificantly smaller than the cylinder used to select isolated 
galaxies and satellites. This is because filaments are not viri- 
alised structures and their intrinsic velocity dispersion about 
the Hubble flow is much lower than that inside a halo (for ex- 
ample, the scatter about the Hubble flow among the galaxies 
surrounding the Local Group is a mere 85 km s _1 , or as low 
as 40 km s -1 if galaxies in side virialised groups are excluded; 
iKarachentsev et al.i r2003), and therefore a much smaller ad- 
ditional velocity is required to account for peculiar velocities 
on top of the Hubble component of 300 km s _1 correspond- 
ing to the radial dimension of the cylinder. We calculate the 
PA and axial ratio of the distribution on the sky of these 
surrounding galaxies by diagonalising the moment of iner- 
tia tensor relative to the primary galaxy. We have used the 
mock catalogues to confirm that this procedure reliably re- 
covers the three-dimensional PA of the LSS surrounding the 
primary (see Appendix |E)) . 

Figure [9] shows the distribution of angles between the 
PA of the great circle connecting the primary and the satel- 
lite and the PA of the LSS surrounding the primary ( "LSS 
angles"; as with the disc angles, these are folded into the 




120 
100 
80 

60 : 

40 

20 



15 30 4 5 60 7 5 9 00 1 5 3 4 5 6 7 5 9 
LSS angle [degrees] LSS angle [degrees] 

Figure 9. (Left) Cumulative distribution of angle between the 
axis of the local LSS and the location of the satellite ("LSS an- 
gle"). The thick black solid/red dotted/green dot-dashed/bluc 
dashed lines refer to the distribution of satellites around all 
primaries, early-type primaries, intcrmcdiate-type primaries and 
late-type primaries respectively. The thin dotted line shows the 
distribution expected if satellites are isotropically distributed. 
(Right) Differential distribution of the LSS angles. The error bars 
are determined by bootstrap resampling of all primary galaxies. 
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Figure 10. Mean LSS angle of satellites as a function of their pro- 
jected separation from the primary. Black plusses, red diamonds, 
green triangles and blue squares refer to satellites of all primaries, 
early-type primaries, intermediate-type primaries and late-type 
primaries respectively. Symbols are offset in radius slightly for 
clarity. 



range 0-90°). We present the results from KS tests and the 
mean LSS angles in Table [5] 

The satellites of each population of galaxies are consis- 
tent with being isotropically distributed with respect to the 
LSS. However, all samples except that around late-type pri- 
maries have mean LSS angles less than 45° . Further data are 
required to determine if this hint of an alignment is real. In 
Figure 1101 we plot the mean LSS angle as a function of ra- 
dial separation from the primary. Although each individual 
point is consistent with isotropy, the satellites at large sepa- 
rations all have mean LSS angles less than 45° (as there are 
by definition no satellites of late-types in this bin, this may 
explain why the late- types also show no hint of LSS align- 
ment). Further information may be gained by plotting the 
LSS alignment as a function of satellite type (Figures llll and 
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Figure 8. Effects of varying the selection parameters on the results and sample size. In each plot, the mean disc angle for 
early/intermediate/late-type primaries is shown as the red diamond/green triangle/blue square symbols with error bars, while the 
sample size for the early- and late-type samples are shown as histograms below and above. The scale for each histogram is shown on the 
right side of each plot. The histograms for the intermediate- type samples show identical trends, (a) The effects of replacing the fiducial 
NED plus photometric redshift ("NED+Pz") method for dealing with violators with either a pure photometric redshift method ("Pz") 
or with a cut on N w - lo \. (b) The effects of changing / sa t- (c) The effects of limiting primaries to be at least m lnneI magnitudes brighter 
than the survey limit and of limiting primaries to be at least projected radii Router from the edge of the photometric survey footprint 
and -R sa t from the edge of the spectroscopic survey footprint. 



I12p . Although the number of early- and intermediate- type 
satellites is too small to draw any conclusions, the late-type 
satellites that are found at large radius are aligned with the 
LSS at the 2a level. 

Given the orientation of the primary galaxy and the 
LSS, we now investigate their relative alignment. The orien- 
tation of a disc galaxy is determined by its angular mo- 
mentum, which originates from tidal torques due to the 
surrounding material. Analytic arguments and cosmological 
simulations suggest that this angular momentum (and there- 
fore the disc spin axis) aligns with the intermediate axis of 
the surrounding mass distribution, and such alignment has 
been measured for d isc galaxies in the supergalactic plane 
(|Navarro et al.ll2004 ) and for galaxies in SPS S and 2dFGRS 
on the surfaces of voids l|Truiillo et al.l [20061 ). The orienta- 
tion of an early-type galaxy is determined by its anisotropic 
velocity ellipsoid, as is that of its halo; therefore, the two 
are expected to be align ed, and preferentially ali gned with 
the large scale structure (jBailin fc Steinmetj2005h . This has 



been measured for B r ightest Cluster Galaxies (BCGs) at low 
jArgvres et al.|[l986l; lLambas et al.lll988l ; iMuriel fc Lambasl 
1 19891 ) and high (|Donoso et al.ll2006l ) redshift, but not for 
field early-types. We directly compared the PA of our pri- 
mary galaxies to that of their local LSS. For this comparison, 
we use all isolated galaxies that pass both the "Disc" and 
"LSS" quality cuts, regardless of whether they host satellite 
galaxies; our results are unchanged if we only include those 
that host satellite galaxies. The distributions are shown in 
Figure [T3l and the associated mean angles and KS test prob- 
abilities of being drawn from an isotropic distribution are 
given in Table [5] There is a detection of alignment between 
the orientation of isolated early-type galaxies and the sur- 
rounding LSS at 99% confidence, a strong alignment for iso- 
lated intermediate-type galaxies at 99.94% confidence. We 
do not detect a significant alignment for isolated late-type 
galaxies. The samples containing all isolated galaxies and 
only those with satellites are consistent with each other. 
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Table 6. Anisotropy of Local Group Satellites 
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Figure 11. As in Figure [9] but separated by the classification of 
the satellite galaxy. 
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Figure 12. As in Figure ITOl but separated by the classification 
of the satellite galaxy. 
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Figure 13. Cumulative (left) and differential (right) distribu- 
tions of angles between the major axes of isolated galaxies and 
their surrounding LSS. The sample contains all isolated galaxies 
that pass the "Disc" and "LSS" sample quality cuts, regardless 
of whether they host any satellite galaxies. Red/dotted lines refer 
to early- type galaxies while blue/dashed lines refer to late- type 
galaxies. 



Parameter 


Milky Way 


M31 


Mean disc angle [°] 


49.1 


54.0 


Median disc angle [°] 


47.7 


53.4 


Polar fraction 


0.57 


0.60 



3.3 Satellites in the Local Group 



Many satellites of the two dominant galaxies in the Local 
Group, the Milky Way and M31, appear to li e on planes 
that are highly inclined to their parent discs l|Metz et al.l 
2007, and references therein). In order to determine how the 
anisotropy of the SDSS galaxies compares to the anisotropy 
around the Local Group spirals, we have determined the sig- 
nal we would have recovered around both the Milky Way and 
M31 if they had fallen into our sampk-0- Based on the dis- 
tribution of primary and satellite magnitudes in our sample 
(Figure Oi), we use all satellites with absolute magnitudes 
within 5 magnitudes of their primary; for the Milky Way this 
consists of the LMC and SMC, and for M31 this consists of 
M33, IC 10, M32, and NGC 205. The 3D locations of these 
satell ites with res p ect to their parent galactic disc are taken 
from iMetz et al.l l|2007t h using the iMcConnachie fc Irwinl 
(|2006h parameters for the M31 satellites. We select 20000 
random viewing directions isotropically distributed about 
each galaxy and calculate the disc angle for each satellite 
from each viewing direction. We calculate the mean disc an- 
gle of the satellites averaged over all viewing angles where 
the projected axis ratio of the parent disc is less than 0.8, 
to provide a direct comparison to the SDSS sample. 

The resulting anisotropies are listed in Table [6] The 
satellites of the Local Group spirals show polar distributions 
with mean disc angles of 49.1° and 54.0° for the satellites 
of the Milky Way and M31 respectively. We do not mea- 
sure such a polar distribution around late-type disc galaxies 
in SDSS. However, we do measure a minor-axis alignment 
of this magnitude among those satellites much fainter than 
their primary (see Figure [5p) and most of the Local Group 
satellites above meet that description. Therefore, perhaps 
the dependence of satellite anisotropy on degree of primary 
dominance may explain the discrepancy between the satel- 
lites of the Local Group spirals and the results from galaxy 
redshift surveys, or perhaps we are simply victims of coin- 
cidence and small number statistics. 



4 COMPARISON WITH PREVIOUS RESULTS: 
THE EFFECTS OF ENVIRONMENT 

Our results indicate that the satellites of isolated early-type 
galaxies show a preference for lying near the major axis of 
the primary, while those of intermediate-type galaxies may 
lie near the disc plane and those of late-type galaxies are 



6 In fact, neither the Milky Way nor M31 would fall into our 
sample if they were observed in a redshift survey because the 
presence of each would violate the isolation criteria around the 
other. Our isolation criteria are required to be so strict in order 
that systems intrinsically less isolated than Local Group galaxies 
are not mistakenly included. 
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isotropically distributed. The result of combining these pop- 
ulations gives a distribution that is purely a function of the 
morphological mixture in the sample of isolated galaxies; in 
our case, the early- and intermediate-type galaxies domi- 
nate over the late-type galaxies and therefore we find that 
satellites tend to exhibit a net major axis alignment. 

It is interesting to compare our results to those of previ- 
ous studies. Our results for late-type galaxies agree wi th pre- 
vious detections of a n isotropic satellite distribution ||SL04| ; 
lAZPKt lAPPZt |Y06| ). and our results for early- type galax- 
ies agree wi t h previous det ections of major-axis alignment 
i|SL04l ; lYOd ; lAPPZt |AB07| ). Our intermediate- type galax- 
ies would have been identified as part of the "early" -type 
sample by those studies that used colour to classify galax- 
ies, and would have been split between apparent ellipticals 
and apparent disc galaxies depending on inc lination by those 
studie s that classified galaxies by eye (see iBailin fc Harris! 
l2008bT ). As the anisotropy around intermediate- type galax- 
ies and early-type galaxies is similar, the former situation 
would not have affected the measured distribution, while in 
the latter case the relatively small number of contaminat- 
ing intermediate-types would not have strongly affected the 
measured distribution around l ate-types. F i nally , our results 
for the full sample agree with |SL04 |B05|. lYOrj and IAB07I . 
who all found that the full sample shows major-axis align- 
ment. 

The main disagreement between our results and previ- 
ous studies are with those studies that fou nd a polar align- 
ment o f satellites around late-type galaxies l|Holmberdl 19691 ; 
IZSFWl ). However, there are physical sub-classes of systems 
for which our results are consistent with a polar alignment 
(i.e. although our results in these regimes are consistent with 
isotropy and we therefore do not claim detection of a polar 
alignment, the mean disc angles are sufficiently larger than 
45° that they are also consistent with a polar alignment): 
satellites at intermediate separation from their primary, and 
those that are much fainter than their primary. The polar 
alignment found by these studies may be explaine d if they 
were dominated by such satellites; indeed, IZSFWl detected 
their polar alignment for satellites with similar projected 
separations as our intermediate separation bin, and domi- 
nance by relatively faint satellites may explain the alignment 
seen around the Milky Way and M31 (see 5 I3.3[I . 

To determine the effects of the different selection crite- 
ria on the measured anisotropy, we adopt the criteria from 
the previous studies (as given in Table [1} to select corre- 
sponding samples of galaxies from SDSS DR6 and to mea- 
sure the disc angle distribution. The mean value of the disc 
angle determined by the previous studies and the value we 
derive using identical selection criteria are shown in Fig- 
ure Q3] as a function of / n0 n-dom, the fraction of satellites 
estimated to lie around non-isolated primaries according to 
the mock catalogue analysis (see Table [3]). For simplicity, in 
this Figure only we separate galaxies into two classes, using 
the location on the CMD. More details about the compari- 
son are presented in Appendix [F] This Figure confirms that 
if we select galaxies according to the criteria used in each 
previous study, we recover their results. 

The angular distribution found using each sample is 
quite similar. However, the selected galaxies lie in very dif- 
ferent environments. Our selection criteria have been fine- 
tuned using mock catalogues to select satellites of isolated 



primaries. While most previous authors have also imple- 
mented selection criteria aimed at identifying isolated pri- 
maries, the results of §[23] indicate that they have had vari- 
able success. In particular, the fraction of satellites around 
non-isolated primaries is over 50% in several of the previous 
studies. The satellites in these systems should be consid- 
ered group members r ather than satellites of iso l ated g alax- 
ies, as should thos e of lYOrJ . TFaltenbacher et all (|2007h . and 
IWang et al.l (|2008l ) . whose selection criteria were tuned to 
find associated galaxies with no constraints on whether the 
largest galaxy in each group is isolated. 

The combination of our results and those of previous 
studies constrain the environmental dependence (or lack 
thereof) of the satellite distribution. Satellites surrounding 
spheroidal galaxies show the same major-axis alignment re- 
gardless of whether that spheroid is isolated or is at the 
centre of a group. Similarly, the satellites surrounding iso- 
lated late-type galaxies are as isotropically distributed as 
the members of groups that surround late-type galaxies. Al- 
though previous studies have not identified intermediate- 
type galaxies, we can select galaxies using their criteria 
and examine the distribution of those satellites around 
intermediate-type galaxies. For example, when using the 
sample generated by the |B05| S2 criteria, which lies at 
the far right side of Figure [14] and contains almost 85% 
group members, the satellites of intermediate-type galax- 
ies have a mean disc angle of 41? 5 ± 0?9, in good agreement 
with the planar alignment tentatively detected around iso- 
lated intermediate-type galaxies using our fiducial criteria. 
Therefore, we conclude that group-specific processes are not 
responsible for the angular distributions of their member 
galaxies, but rather processes that also apply to the satel- 
lites of isolated galaxies. 



5 SUMMARY AND DISCUSSION 

We summarise our results as follows: 

• Satellites of isolated early-type spheroidal galaxies lie 
preferentially along the major axis of the galaxy. The degree 
of alignment increases slightly for satellites that are brighter 
relative to their primary. 

• Satellites of isolated disc galaxies appear to have dif- 
ferent angular distributions depending on the colour of the 
disc. Satellites of red discs (intermediate-type galaxies) show 
hints of lying preferentially near the disc plane, with intrin- 
sically or relatively fainter satellites showing a tendancy to 
stronger anisotropy. Satellites of blue discs (late-type galax- 
ies) are distributed isotropically. 

• Late-type satellites that are found far from their pri- 
mary show preferential alignment with the surrounding LSS 
(i.e. filaments) (Figure [12]). 

• Isolated early- and intermediate-type galaxies show an 
alignment with the surrounding LSS. This alignment is 
strongest for intermediate-type galaxies (with a KS test sig- 
nificance of 99.94%). 

• The angular distribution of group members about the 
BGG is very similar to the distribution of satellites around 
an isolated galaxy of the same type. 

• Great care must be taken in order to select truly iso- 
lated galaxies and their satellites in galaxy redshift sur- 
veys. Unless the region immediately surrounding the pri- 
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mary is devoid of galaxies too large to be considered satel- 
lites (whether or not those galaxies have been observed spec- 
troscopically) , the sample will be dominated by group mem- 
bers rather than isolated galaxies. 

These results provide us with the foundations on which 
we can build our understanding of the mass distribution 
in and around galaxies and raise a number of interesting 
questions. What can we learn about the role of dynamical 
effects in driving preferential alignments of satellites? What 
is the nature of these effects and does it depend on the 
morphological type and history of the galaxy? What role is 
played by the host dark matter halo? What role is played 
by the larger scale environment? Our results allow us to 
begin to address these questions. 

The most straightforward interpretation of the prefer- 
ential alignment of satellites at large projected radii with the 
surrounding LSS is that this is a signature of the anisotropic 
infall of satellite galaxies along filaments. This interpretation 
is favoured by the following observations: 

• the mean angle between the projected positions of satel- 
lites and the surrounding LSS (i.e. the LSS angle) is similar 
around different types of galaxies; 

• the alignment is present for only the more recently ac- 
creted satellites, i.e. those that are most distant and those 



blue late-type satellites for which halo-specific transforma- 
tion mechanisms have not yet had time to operate; and 

• the alignment is distinctive when compared to the align- 
ment of the satellites with respect to their primary. 

This evidence argues in favour of an origin that is imposed 
by the larger scale environment rather than one driven 
by the primary galaxy. In other words, it is improbable 
that the dynamics of the most recently accreted satellites 
will be significantly affected by any process internal to 
the primary's dark matter halo. This is in agreement with 
the results of cosmological simulations l|Knebe et al.l |2004| ; 
IZentner et al. l l2005l : lLibeskind et alj|2005h . 

The relationship between the orientation of a galaxy 
and its surrounding LSS can be understood in terms of the 
relationship between the orientation of the galaxy's dark 
matter halo and its surrounding LSS and the galaxy's ori- 
entation within its dark matter halo. Cosmological iV-body 
simulations predict that dark matter haloes a re strongly tri- 
atrial systems in the absence of baryons (e.g. lAllgood et al.l 
2006), and these haloes tend to align with their major axes 
along the large scale filaments and their minor axes perpen- 
dicular to filaments (jBailin fc Steinmetdl2005h . 

However, the presence of baryons can have a dramatic 
effect on the shapes and internal alignments of dark mat- 
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ter haloes, with interesting consequences for the haloes of 
disc galaxies. A number of studies have shown that cool- 
ing baryons at the centre of a dark matter halo tend to 
circularise the orbits of dark mat ter particles, modifying 
the halo's inner mass profile (e.g. iGustafsson et all 120061 ) 
and reducing the ellipticity of the halo's isode n sity surfaces 
(e.g. iDubinskil Il994l ; iKazantzidis et all 120041 ). iBailin et ail 
(2005) examined the structure of the host haloes of several 
disc galaxies that formed in high resolution cosmological 
iV-body hydrodynamical simulations, and discovered that 
haloes consisted of two distinct regions. The inner halo is 
flattened along the disk axis, while the orientation of the 
outer halo is unrelated to that of the inner halo and is un- 
affected by the presence of t he luminous galaxy. 

The lBailin et~ail l|2005l ) result is interesting because it 
implies that the major axis of the inner halo around a typ- 
ical disc galaxy is aligned with the major axis of the light 
distribution, while the major axis in the outer region is inde- 
pendent of the light but is aligned with the LSS. If satellite 
galaxies are more common along the halo major axis, then 
this "twisting" of the halo should be evident in the distri- 
bution of satellites around disc galaxies. Indeed, we observe 
that it is the outermost satellites around all types of galax- 
ies that are preferentially aligned with the LSS, while it is 
the inner satellites around intermediate-type disc galaxies 
(which correspond most closely to the relatively red c oncen- 
trated simulated discs studied bv lBailin et al. (2005)) that 
are found preferentially in the disc plane. However, the sig- 
nificant relative alignment between intermediate-type galax- 
ies and the LSS argues that some residual halo-LSS align- 
ment remains. 

In the case of early-type galaxies, the orientation of 
the galaxy with respect to its dark matter halo has not 
been studied explicitly in a cosmological context. However, 
we would expect that the shapes of both their stellar and 
dark matter components are supported by their anisotropic 
velocity ellipsoid. Consequently, we would expect that both 
the galaxy and its dark matter halo will share the same 
orientation, and therefore the galaxy will tend to align 
with the surrounding LSS (IBailin fc Steinmet j [2005). as is 
observed. 

How should we interpret the preferential alignments 
of satellites around isolated galaxies? Can we determine 
whether the alignments are imprinted by the dynamical ef- 
fects lofthegala^rorthehostd^xk matter halo? 

_ Agustsson fc Brainerdl l|2006l ). IAB07I and iKang et al.l 
studied the angular distribution of satellite galaxies 
in cosmologica l sim ulati ons selected according to the criteria 
of |B05| . IABOTI and lYOfj respectively (note that all of these 
criteria select samples dominated by groups). The orienta- 
tion of a mock galaxy must be assumed and so these authors 
explored different assumptions about how primary galaxies 
are oriented with respect to their dark matter haloes and 
the larger scale environment. They found that if the pri- 
mary galaxy is a spheroid whose principal axes are perfectly 
aligned with those of its dark matter halo, then satellites in 
these systems tend to show a major-axis anisotropy that is 
stronger than observed in groups whose BGG is an early- 
type galaxy. However, if there is a small offset between the 
principal axes of the galaxy and its halo, as may arise if the 
galaxy aligns with the halo's angular momentum rather than 



its minor axis, then the anisotropy is of the same order that 
is observed. Because the alignment of satellites relative to 
isolated early-types is identical to that seen in galaxy groups 
whose BGG is an early-type, we may therefore conclude that 
the principal axes of isolated early-types are similarly well 
aligned with those of their dark matter haloes (~ 20°). As 
we argued above, this is in accord with our expectation that 
the dynamics of baryons and dark matter in collisionless 
ellipsoidal systems are similar. 

The theoretical situation around discs is less clear. If 
the disc is oriented perpendicular to its halo's angular mo- 
mentum, it is simply a special case of an oblate spheroid, 
and therefore these studies predict that its satellites will ex- 
hibit major-axis alignment. This is consistent with what we 
observe around red discs, but in stark contrast to what is 
seen around blue discs. If, on the other hand, the angular 
momenta of galaxy discs align with the intermediate axis 
of the su rrounding mass distr ibution, as seen in the simu- 
lations of iNavarro et al] (|2004l ). then satellites in these sys- 
tems show no preferential alignments, as we observe around 
blue discs. If we interpret disc colour as a measure of how 
long the baryonic material has been part of the luminous 
galaxy, then perhaps red discs have had more time to come 
to equilibrium with their halo, while most of the material 
in blue discs has been acquired more recently and retains a 
memory of the external tidal torques. 

However, these explanations are still largely speculative. 
A theoretical analysis that takes into account the detailed 
dependences of the satellite distribution, the alignment of 
satellites with the LSS, and the differences between early- 
type spheroidal galaxies, intermediate- type red disc galaxies, 
and late-type blue disc galaxies must be performed to de- 
termine whether the distribution of satellites is determined 
predominantly by the orientation of the halo or if dynamical 
processes within the halo are important. We are in the pro- 
cess of performing such an analysis (Power et al., in prepa- 
ration). 
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mass, and proposed a parameterisation for the variation of 
the average total mass-to-light ratio with halo mass, 



/Mr 11 / M \ 



(Al) 



Here the free parameters correspond to Mi, the character- 
istic mass for which the mass-to-light ratio in bj is equal 
to (M/L)o, and 71 and 72 which determine the behaviour 
at the low- and high- mass ends of the of the m ass function 
respectively. We follow lYang et all (|2004T ) (|Y04T l in adopting 
M 1 =10 w - 94 h- 1 M Q , (M/L) =124 h (M/L) in bj, 71 =2.02, 
and 72 =0.30. 

The characteristic luminosity L* is parameterised in a 
similar manner; 



M 1 (M\ ,, .\fM\ ~ 71 / M 



(A2) 



where M2 is a characteristic mass and 73 determines the be- 
haviour at the high-mass end of the mass function; a follows 



a = au + ?7log(Mi5) 



(A3) 



where M15 is the mass of the halo in units of 10 1 5 fo 1 Mq 
and Q15 and 77 are free parameters. We follow Iy53 in adopt- 
ing A/ 2 =lO 12 O4 /i- 1 M , 73 =0.72, r/=-0.22, and a 15 =-l.l. 

Expression ! All and IA2l allow an expression for (L) / (M) 
to be derived, from which $* is deduced; 
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(M) 



<M) 
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$*(M) = 
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HL\M) — dL = $*— F(q + 2) (A4) 
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Here r(i) and T(x,a) are the Gamma and Incomplete 
Gamma functions respectively; formally these are expressed 



F(x) = / t*' 1 exp(-t)dt 
Jo 

and 



r(x,a) 



t x exp(—t)dt 



(A6) 
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We note that IyoI denote the average mass-to-light ratio by 
(M/L) (see their equation 2). We prefer (M) / (L) because 
the meaning is clear - the average luminosity associated with 
a halo of mass M is (L) and so the average mass-to-light 
ratio is (M) / (L). If we adopt (M/L), this means that 



(t><">- 



M M 
$(L\M)-^dL = $*4^r(a). (A8) 
L £,* 



APPENDIX A: THE CONDITIONAL 
LUMINOSITY FUNCTION : PARAMETERS 



lYang et all (120031 ) deduced a functional form for the vari- 
ation of the mass-to- light ratio with dark matter halo 
mass by comparing the lSheth fc Tormenl (|l999l ) dark matter 
halo mass funct ion with the Schechter luminosity function 
l|Schechterlll976l ). They noted that the mass-to- light ratio 
must increase (decrease) with decreasing (increasing) halo 



This produces an expression for $* that is quite different 
from equation IA5I and which does not recover the correct 
behaviour of quantities such as (TV) (M). 

Having deduced the form of $*, we can compute the 
"conditional luminosity function" $(L|M), 



$(L\M)dL = i- (^j-^j exp{-L/L*)dL 
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The upper left-hand panels of Figure IA1I show how L* 
and <E>* vary with halo mass, while the right-hand panel 
shows the variati on of <&(L|M) with luminosity at a fixed 
halo mass for the choice of 2dFGRS parameters. 



"hard" lower mass cut-off of M m j n = 10 9 h~ 1 M t r 1 below 
which haloes cannot host galaxies - galaxy formation is 
suppressed in these haloes following cosmological reionisa- 
tion. 



APPENDIX B: POPULATING DARK MATTER 
HALOES WITH GALAXIES : DETAILS 

We perform a suite of cosmological iV-body simulations 
and constructed catalogues of dark matter haloes at z~0. 
Dark matter haloes are identified using a friends-of-friends 
(FOF) algorithm with a linking length of b = 0.2 times 
the mean interparticle separation. For each of the groups 
identified in this way we compute the virial mass Migo, 
defined as the mass of the spherical overdensity that is 180 
times the critical density of the Universe at z=0. In the 
following discussion, we define a halo's mass M to be its 
virial mass Miso rather than Mfq f , the mass of the FOF 
group; this is required by the IyoI prescription. 

The minimum halo mass M m i n that is "reliably" recov- 
ered in each of the simulations governs the minimum lu- 
minosity, I/min, that is used in constructing the mock cat- 
alogues. L m i n defines the threshold luminosity fainter than 
which there are no galaxies. M m i n is the halo mass above 
which we expect the mass function to be unaffected by 
finite numerical resolution; below this threshold the num- 
ber density of haloes tends to be suppressed relative to the 
number density they would have in the limit of infinite nu- 
merical resolution. Previous studies have examined how the 
mass function is affected by finite mass and force resolu- 
tion, time-stepping accuracy and starting redshift, as well 
as the influence of the grou p-finding algorithm used to iden- 
tify d a rk matter halo e s (e.g . Ijenkins et al.|[200ll ; iLukic et al.l 
120071 ). I Jenkins" et alj (|200ll ) performed careful convergence 
tests and found that mass functions constructed from FOF 
groups are adversely affected by numerical effects below a 
halo mass equivalent to ~ 20 particles. In this work we adopt 
a more conservative lower mass limit of 50 particles, to en- 
sure that the mass function of halo es in higher density re- 
gions is converg ed (see Chapter 3 of |Power|[2003T ): this gives 
M min =50m par t in Table [2] 

Having determined M m i n , we estimate L m m using the 
"conditional pro babili ty distribution" P(M\L) (see right 
panel of Figure 1. |Y04) . Z/ m i n is a critical parameter because 
it fixes the halo occupation number, the average number of 
galaxies per halo of mass M; 

(N) (M) = $*r(Q + 1, L min /L*). (Bl) 

The variation of (N) (M) with halo mass for the lY04l 2dF- 
GRS parameters and L m i n — 1.1 x 10 9 /i~ 2 L-0 is shown in 
the bottom left-hand panel of Figure I All This L m i n is ap- 
propriate for the M m j n in the mock catalogues A to E. 

Note the important role played by the ratio L m i n /L* 
in equation IB1I which controls the number of galaxies per 
halo. At fixed L m i n it increases dramatically as halo mass 
decreases, leading to low mass haloes containing one "cen- 
tral" galaxy at most, and as L m i n decreases, the number of 
galaxies per halo increases. The number of galaxies per halo 
of mass M is Poisson dis tributed wit h a me an of (N) (M). 

We note also that lYang et all |2003h introduced a 



Having determined the number of galaxi es ho sted by a 
halo, we must assign luminosities. We follow |Y04| and give 
special status to the central galaxy by assuming that it is 
the brightest in the halo, with an average luminosity 



(L c ) = $*Z*F(Q + 2,ii/L* 



(B2) 



The luminosity L\ is a function of halo mass and is chosen 
such that 



l>*r(a + l,Li/L*) 



(B3) 



when choosing the central galaxy luminosity, we assume 
that L c is a random variable drawn from $(L|M) for the 
range of luminosities L > L\. The remaining N — 1 galaxies 
within the halo are assigned luminosities in the range 
Lmin < L < Li, drawn at random fro m th e luminosity 
function (the "intermediate" approach of |Y04 ). 

The penultimate step involves assigning morphological 
types to each mock galaxy; this is done by defining a func- 
tion /i a t e (L,M) that specifies the fraction of galaxies with 
luminosity L in haloes of mass M that are late-type. This 
function can be expressed as the product of functions, 



/iate(i, M) — g(L)h(M)q(L, M), 



where 
q(L,M) = 



1 if g{L) h(M) < 1 

WTHM) X9(L)h(M)>l 



g(L) = 



and 



&i.f(£) f~*(L\M)n(M)dM 
-I(L) J °° ${L\M)h(M)n{M)AM' 



Iog(M/M„) 



(B4) 

(B5) 
(B6) 

(B7) 



h(M) = max 0, min , . , , , , , , 

1 \\og{M b /M a ) 

Here n(M) is the halo mass function (|Sheth fc Tormenl 
1999); $i a tc(^) and $(i) correspond to the observed lumi- 
nosity functions of the late-type and entire galaxy samples 
respectively; and M a and Mb are free parameters defined 
as the masses at which h(M) takes on th e values and 
1 respectively. Ivan den Bosch et al l (|2003l ) demonstrated 
that this parameterisation allowed the galaxy population to 
be split into early- and late-types such that the respective 
luminosity functions and clustering properties could be 
recovered. We follow |Y04| in adopting Af a = lO 17 - 26 ^" 1 ^^ 
and Mb = lO lo ' 86 /i -1 M0. Formally we assign morphological 
type by drawing a random number R that is uniformly 
distributed between [0,1] and comparing it to f\ ate (L, M). 
If R < /iate(i,Af), the galaxy is designated late-type, 
otherwise it is early-type. 

The final step involves assigning phase space coordi- 
nates (i.e. positions and velocities) to each of the N galaxies 
within the halo. The brightest central galaxy is associated 
with the most bound particle of the halo and is assigned 
its position and velocity. The remaining N-l galaxies can 
be treated in a variety of ways. For the purposes of this 
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Figure Al. (LEFT) Variation of the conditional luminosity function par amete rs L* (upper panel) and <&* (middle panel), and the halo 
occupation number (N) (M) (bottom panel) with halo mass. We adopt the lY04l choice of 2dFGRS parameters and a minimum lumin osity 
^min = 1-1 x 10 9 /i — 2 Lq. (RIGHT) Variation of the conditional luminosity function with luminosity at fixed halo masses, for the I Y04l 
parameters. 



study, in which our main concer n is t esting; the reliability of 
our selection criteria, we follow IyoI in randomly sampling 
dark matter particles from the FOF group (their "FOF ap- 
proach"). More sophisticated approaches, in which we ex- 
plicitly track the merging history of individual haloes, will 
be essential for future work, especially with regards to kine- 
matics (Power et al., in preparation). 



APPENDIX C: GALAXY CLASSIFICATION 

Our pri mary classification method is that of lBailin fc Harris] 
l|2008bl ). which has been validated using high-q uality imag- 
ing f rom the Millennium Galaxy Catalogue (Lis ke et alj 
2003). However, given the qualitative difference between our 
results for the different subpopulations, it is important to 
confirm that the difference seen is not an artefact of the 
galaxy classification scheme. We examine how the anisotropy 
of the satellite distribution varies using other methods: 

(i) Inclination-corrected location on the colour- 
magnitude diagram (C MD F ), which is strongly bimodal 
( Bail in fc Harris] l2008bh . Galaxies are considered "early" if 
they are redder than CMD F = —0.05, and "late" if they are 
bluer. 

(ii) Spectroscopic Principal Component Analysis (PCA) 
eClass parameter. Galaxies are considered "early" if they 
have eClass < —0.07, otherwise they are considered "late". 

(iii) The inclination-corrected global concentration of 
the light profile, Cnorm. The distribution of ga laxy con- 
centrations is trimodal (jBailin fc Harris! l2008al ). We la- 
bel the "Elliptical" (hi gh-C n orm and high-b/a) region from 
iBailin fc Harris! |2008al ) as "early" , their "Disk" (low-C nor m) 
region as "late" , and all other galaxies as "intermediate" . 

It should be noted that these measurements are completely 
independent: the CMD location is based on global photom- 



etry, the PCA analysis is based on spectroscopy, and the 
concentration is based on the distribution of the light pro- 
file. 

The results using these alternative classification 
schemes are shown in Figure [Cll The mean disc angle and 
KS test probability that each sample is drawn from an 
isotropic distribution are also given. The major-axis distri- 
bution around early-type galaxies and a distribution con- 
sistent with isotropy around late-type galaxies is seen us- 
ing every method. Intermediate-type galaxies are red with 
intermediate concentrations; the galaxies with red CMD F 
and with intermediate Cnorm show the same major-axis 
distribution as around the intermediate-type galaxies of 
IBailin fc H arris ( 2008i3) ■ Therefore, although the magnitude 
of the measured anisotropy varies at a ~ la level, the de- 
tected anisotropy cannot be simply a galaxy classification 
artefact: satellites of early-type and late-type galaxies have 
different angular distributions. 

We note that the classification of some of our galaxies is 
uncertain. Because the galaxies that constitute our primary 
sample are typically more luminous and more isolated than 
typical SDSS spectroscopic galaxies, they provide a biased 
sample of parameter space. In particular, several of our pri- 
mary galaxies have Cnorm < 1 (i.e. they have low concentra- 
tions) but have CMD F > -0.05 (i.e. they are red): 4.3% of 
all primaries and 12.0% of primaries classified as Late-type 
fall into this region of parameter space, c ompared to just 
1.8% o f the visually-classified galaxies in IBailin fc Harris! 
(|2008bl ) and 3.6% of those classified as Late. Given that 
the anisotropy of the satellite distribution differs between 
galaxy classes, and shows the strongest difference between 
the Late and Intermediate types, examining the anisotropy 
around the galaxies in this region of parameter space can 
provide insight into their nature. 

If we separate our Late types into red and blue sub- 
classes (divided at CMD F = —0.05, as above), we find that 
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the mean disc angle around the blue subclass is 46.7 ± 2.2°, 
consistent with isotropy and with the results from the full 
Late-type sample. However, the mean disc angle around the 
red subclass is 37?3 +5 i°, exhibiting major-axis alignment 

—5. 5 

consistent with the results from the Intermediate-type sam- 
ple. This suggests that the red low-concentration galaxies 
more properly belong in the Intermediate classification. 



APPENDIX D: ANGULAR DISTRIBUTION OF 
SATELLITES IN THE MOCK CATALOGUES 

To confirm that our measurement of an anisotropic distri- 
bution of satellite galaxies is due to an intrinsic anisotropy 
rather than an artefact of our method, we have performed 
an identical analysis on the mock catalogues, whose satellite 
distributions are, by construction, isotropic. We use mock 
catalogues generated from the five independent 100 h~ x Mpc 
simulations in order to account for cosmic variance. The dis- 
tributions are plotted in Figure IDll and the statistical mea- 
sures of anisotropy are listed in Table IDll 

The level of anisotropy that we measure from 
isotropically-distributed satellites in the mock catalogues is 
small. Even when the isotropic KS test probabilities in the 
mock samples are low, the deviations are not systematic: the 
mean disc angle almost always deviates from 45° by less than 
2° the polar fraction never deviates from 0.5 by more than 
4%. These are much smaller than the anisotropies that we 
detect around early-type galaxies in the observational sam- 
ple, confirming that our detection of anisotropy cannot be 
explained by intrinsically isotropically distributed satellites. 



APPENDIX E: DETERMINATION OF THE 
LARGE SCALE STRUCTURE AXIS 

The physical environment of a galaxy is best described by 
the region in which the presence of the galaxy predicts the 
presence of other matter; the radial extent of this region is 
characterised by the correlation length ro. For the global 
~ L* galaxy population, ro ~ 4-6 h~ x Mpc jNorberg et al.l 
2002). However, the isolated galaxies that constitute our 
sample are, by construction, much less clustered than aver- 
age; for example, Hl-selected galaxies, which are much less 
likely than average to have large ne arby neighbours , have 
a much smaller r p ~ 3.3 h~ r Mpc (|Basilakos et al.l 120071 ; 
iMever et al1l2007l ). For our very isolated sample, 3 h 1 Mpc 
is a reasonable radius in which to characterise the large-scale 
environment of each galaxy. 

We therefore determine the PA of the large scale struc- 
ture around each primary galaxy by diagonalising the mo- 
ment of inertia tensor of the projected positions of all spec- 
troscopic galaxies with projected radii of between 1000 and 
3000ft _1 kpc (thereby explicitly ensuring that there is no 
overlap between the galaxies used to determine the orienta- 
tion of the LSS and those used to evaluate the isolatedness 
of the primary or the satellites themselves) and with veloc- 
ities that differ from that of the primary by no more than 
400 km s _1 (this is larger than the 300 km s -1 Hubble flow 
component in order to account for the peculiar velocities of 
galaxies, which have a dispersion of 85 km s~ in the Local 
Volume; iKarachentsev et al"1l2003l ) . 
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Figure El. Cumulative distribution of the difference between 
the PA of the LSS measured around isolated mock galaxies using 
the "observed" 2D galaxy positions and redshifts versus that mea- 
sured using the known three-dimensional positions of surrounding 
haloes. The solid line indicates galaxies for which both the 2D and 
3D PA is well-defined, while the dashed line also includes haloes 
for which only the 2D PA is well-defined. 

We have used the mock catalogues to confirm that this 
procedure reliably recovers the three-dimensional PA of the 
LSS surrounding the primary. We have taken the known 
three-dimensional positions of all haloes within a spherical 
volume of 3000 h~ x kpc around the halo of each primary 
galaxy in the mock catalogues, constructed and diagonal- 
ized their inertia tensor, and projected the major axis onto 
the plane of the sky. The PA of this axis is then compared to 
the two-dimensional PA inferred from the "observed" galax- 
ies in the mock catalogues. The cumulative distribution of 
the misalignment between the PA determined using 3D posi- 
tions and that inferred from the 2D observables is plotted in 
Figure lETI To improve the statistics, we have included in this 
plot all galaxies that match the isolation criteria, regardless 
of whether they have satellites; however, the results are con- 
sistent if we restrict the sample to just those with satellites. 
The solid line indicates the relative alignment when both the 
3D LSS PA and the inferred 2D are well-defined (i.e. that 
contain galaxies within the defining sphere or cylinder), and 
has a median misalignment of 13.8°. A small fraction of 
galaxies have well-defined 2D LSS PAs, and would therefore 
be included in the observational analysis, but no well-defined 
3D LSS axis because none of the galaxies that lie within the 
redshift-space cylinder lie within the 3D sphere. We account 
for these cases by assuming that their intrinsic 3D LSS PAs 
are isotropically distributed and indicate the alignment of 
the full sample including them as the dashed line in Fig- 
ure Hi] Half of the LSS PAs are aligned to within 15.8°; 
therefore, this procedure successfully recovers the PA of the 
LSS. 



APPENDIX F: COMPARISONS USING 
DIFFERENT SELECTION CRITERIA 

In Figure 1141 we have compared the mean disc angle de- 
termined by previous studies of satellite anisotropy to the 
values we derive using identical selection criteria. 

No previous study has identified intermediate-type 
galaxies as a separate class; we therefore adopt a simple 
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Figure CI. Cumulative (top panels) and differential (bottom panels) distributions of satellite disc angles of primary galaxies classified 
using the following schemes (left to right): CMD F , spectroscopic eClass parameter, and inclination-corrected Petrosian concentration 
parameter C norm . Colours/line styles are as in Figure [4] Mean disc angles and KS test probabilities that the samples are drawn from a 
uniform distribution are given in the top-left (bottom-right) corners of the cumulative plots for the early-type (late-type) samples. The 
statistics for the intermediate- type samples are given below the early-type statistics for C nc ,rm- 



Table Dl. Anisotropy of Satellites in Mock Catalogues 



Parameter Mock A Mock B Mock C Mock D Mock E 
Full sample: 

KS probability 0.28 0.41 0.95 0.94 0.65 

Mean disc angle [°] A4.7+\\ 46.8 ± 1.5 45.3 ± 1.5 44.91^4 45.9+.*;| 

Median disc angle [°] 44.7+2 47A -l'l 44.4+^ AbA+ 2 2 \ 44.1+/^° 

Polar fraction 0.49 ± 0.03 0.53 ± 0.03 0.49 ± 0.03 0.50 ± 0.03 0.50 ± 0.03 
Early-type primaries: 

KS probability 0.30 0.33 0.76 0.46 0.64 

Mean disc angle [°] 46.1+.^° 48.4 ± 2.2 45.0+j g 45.2 ± 1.7 46.2+J® 

Median disc angle [°] 44. 48.8jl5;f 44 - 8 ii'.l 44.0±f'j 

Polar fraction 0.50+° o4 0.55 ± 0.04 0.49 ± 0.04 0.48 ± 0.03 0.50 ± 0.03 
Late-type primaries: 

KS probability 0.10 0.83 0.45 0.40 0.80 

Mean disc angle [°] 43.1+^° 44.9^j 45.7+2 2 ^-^l 45.2 ± 2.2 

Median disc angle [°] 43.7±2 ;* 45.2^;| 43.7±3.2 47.9+J; 4 43.4+^'S 

Polar fraction 0.48 ± 0.04 0.51 ± 0.04 0.48+_£; 4 0.54 ± 0.05 0.48 ± 0.05 



binary classification based on the location of the galaxy on 
the CMD: a galaxy is considered Early-type if its 01 (g — r) 
colour is redder than 

0,1 (ff ~r) = 0.70 - 0.0325(M r - 5 log h + 19) (Fl) 

and Late-type if it is bluer than this threshold. The results 
of Appendix[C]indicate that different methods of classifying 
galaxies may introduce ~ la differences in the measured 
anisotropy. 

We split e ach sam ple by the galaxy type of its pri- 
mar y alt hough IZSFWl and IAZPKI only studied l ate- t ypes, 
and iBOrj did not separate the sample by type. IB05I sug- 
gested that her samples were dominated by systems with 
late-type primaries; in contrast, we find that 46%/58%/37% 



of the primaries we select using her criter i a are cl assified as 
early types according to iBailin fc Harris! (|2008bt ). contain- 
ing 52%/64%/43% of the satellites for samples 1, 2, and 3 
respectively. Using the location on the CMD, the fraction of 
early- type primaries is even higher: 85%/90%/83% contain- 
ing 89%/9 4%/86% of the satellites. 

ISL04I do not quote a mean disc angle; rather, they fit 
the distribution of disc angles 9 to the form 

f(9) = Acos(28) + B (F2) 

and quote the values of A and errors a a- To enable a more 
direct comparison with other studies, we calculate the mean 
disc angle of the associated distribution as 
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Figure Dl. Cumulative (top panels) and differential (bottom panels) distributions of satellite disc angles in the five mock catalogues. 
Colours/line styles are as in FigurefJ] Mean disc angles and KS test probabilities that the samples are drawn from a uniform distribution 
are given in the top-left (bottom-right) corners of the cumulative plots for the early-type (late-type) samples. 



= ( F3 ) 
with uncertainty 

*w = ^, (F4) 

in radians. lAPPZl also do not quote a mean disc angle. We 
have derived the mean and the error of the distributions 
from their plotted histograms, assuming that all satellites 
lie at the central value of the bin the y fall in . The errors are 
calculated by bootstrap resampling. ZSFW do not quote a 
mean disc angle, but it can be derived from the data in 
Table 2 of lZaritskv et al.1 (|l997al ). The errors are calculated 
by bootstrap resampling. 

The only cases where our results deviate from the pre- 
vi ous resu lts by more than 2a are the early-type samples 
of lAPPZl Even in these cases, the sense of the observed 
anisotropy is the same, only the magnitude is different. 
These are cases where there is no mean disc angle quoted 
by the authors, and therefore we have used indirect methods 
to determine the appropriate mean; they also use different 
methods to classify early-type galaxies. We conclude that 
the numerical differences in these cases are not significant. 



